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ABSTRACT 


The  flow  characteristics  of  emulsions,  in  a  porous 
medium,  were  examined  by  studying  the  flow  of  two  low 
concentration  (  r  =  0.15  and  0.30  )  oil-in-water  emulsions, 
through  one  consolidated  core  and  cores  made  from  three 
different  unconsolidated  sands. 

The  results  showed  that  an  emulsion  effectively 
reduces  the  permeability  of  a  core,  when  the  emulsion 
droplets  are  larger  than  the  pore  throats  in  the  porous 
medium.  When  emulsions  plug  a  porous  medium,  the 
relationship  between  flow  rate  and  pressure  drop  is  not 
linear  and  therefore  emulsion  flow  in  a  porous  medium  is  a 
non-Darcy  flow  process.  The  shape  of  the  fractional 
permeability  versus  flow  rate  or  pressure  drop  curve  is  a 
function  of  the  pore  throat  size  distribution  in  the  porous 
medium.  There  is  no  simple  relationship  between  the 
permeability  reduction  and  the  porosity  and  permeability 
of  the  core . 

The  shape  of  the  produced  emulsion  concentration 
profile  was  examined  for  the  0.15  concentration  emulsion. 

The  droplet  size  distribution  of  the  produced  emulsion  was 
examined  in  one  experiment.  The  results  showed  that  the 
degree  of  retention  was  higher  for  the  larger  oil  droplets. 

A  model  based  upon  the  assumption  that  the  droplets 
can  plug  the  porous  medium  and  this  rate  of  plugging  can  be 
modelled  by  a  rate  equation  is  proposed.  The  model  has  to  be 
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compared  with  the  experimental  data  to  evaluate  the  various 
unknown  parameters.  Three  models  were  proposed  in  this 
study  and  one  model  effectively  fitted  the  experimental 
data . 
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CHAPTER  1 


INTRODUCTION 

1.1  BACKGROUND 

In  the  oil  industry,  it  has  been  known  for 
many  years,  from  as  early  as  the  1942  patent  of  Subkow  (l), 
that  emulsion  formation  within  the  reservoir  can  improve  oil 
recoveries.  Since  that  time,  a  number  of  researchers  (2)  - 

(6)  have  reported  a  significant  increase  in  oil  displacement 
efficiency  as  a  result  of  in  situ  emulsification  during 
water  and  caustic  floods . 

Besides  the  possibility  of  naturally  forming  an 
emulsion  within  the  formation  matrix,  a  number  of  processes 
directed  towards  increasing  oil  recovery  efficiencies  by 
injecting  externally  produced  emulsions,  have  been  patented 

(7)  -  (9)»  McAuliffe  (10)  has  reported  the  results  of  the 
only  emulsion  flood  field  test  conducted.  The  wells  in  this 
test  showed  improved  oil  recoveries  and  lower  water-oil 
ratios,  than  the  wells  in  the  surrounding  conventional 
waterflooding  patterns. 

Initially,  interest  in  emulsions  was  primarily  related 
to  obtaining  a  clearer  understanding  of  caustic  flooding. 
This  process  was  recently  reviewed  by  Johnson  (11).  He 
discussed  the  four  principal  mechanisms  for  caustic 
flooding: 

1)  Emulsification  and  Entrapment 

2)  Emulsification  and  Entrainment 

3)  Wettability  Reversal  (  water-wet  to  oil-wet  ) 
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4)  Wettability  Reversal  (  oil-wet  to  water-wet  ) 
of  which  three  involve  the  formation  and  flow  of  emulsions. 

More  recently,  comparisons  and  investigations  on 
caustic,  acid,  water  and  emulsion  floods  of  viscous  oil 
formations  has  received  some  attention  (12)  -  (14).  The 
work  of  Farouq  Ali  et  al  (14)  showed  that  emulsion  slugs 
were  more  effective  in  recovering  oil  than  caustic  or  acid 
slugs,  or  continuous  water  injection.  D'Elia  and  Ferrer  (13) 
and  McAuliffe  (12)  both  reported  improved  recoveries  of 
emulsion  floods  over  waterfloods. 

Emulsion  and  caustic  flooding  methods  have  proven  to 
be  more  effective  than  conventional  non-thermal  methods  of 
recovery  for  moderately  viscous  oils,  because  they  can  in 
part  overcome  the  three  main  factors  responsible  for  low 
recoveries,  which  are: 

1)  the  unfavourable  mobility  ratio, 

2)  the  adverse  surface  and  interfacial  forces  between 
the  oil,  water  and  formation  matrix  and 

3)  the  poor  sweep  efficiency  obtained  during  flooding. 

Principally  emulsion  flooding,  as  an  oil  recovery 

process,  has  an  application  in  only  certain  types  of  heavy 
oil  reservoirs.  The  prime  candidates  are  thin,  relatively 
deep,  high  permeability  formations,  containing  moderately 
viscous  oils.  In  such  situations,  thermal  recovery  may  not 
be  feasible  in  view  of  depth  and  small  thicknesses.  Another 
candidate  is  one  where  the  oil  formation  is  underlain  by 
a  water  bearing  sand  (14). 
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Although  extensive  work  has  been  carried  out  on  the 
oil  recovery  efficiency  of  caustic  and  emulsion  floods, 
relatively  little  has  been  published  on  specific  aspects  of 
emulsion  flow  in  a  porous  medium.  Thus  if  the  processes  of 
emulsion  and  caustic  flooding  are  to  be  carried  out 
economically,  the  properties  of  the  emulsion  and  the  porous 
medium,  and  more  importantly,  the  effects  of  these 
properties  on  the  flow  of  emulsions  through  a  porous  medium 
must  be  clearly  defined. 

With  a  better  understanding  of  emulsion  flow, 
justifiable  engineering  decisions  could  be  made  to  improve 
oil  production.  A  mathematical  model  for  emulsion  flow 
could  also  be  incorporated  into  reservoir  simulators  to 
obtain  more  reliable  estimates  of  such  processes  as  caustic 
and  steam  flooding.  The  end  result  of  the  improved  knowledge 
should  be  to  improve  the  oil  recovery  efficiency. 


1.2  PURPOSE  OF  THE  STUDY 

The  study  of  the  flow  of  macroemulsions  through  a 
porous  medium  was  the  broad  objective.  Specifically  the 
objectives  were: 

1)  to  demonstrate  that  emulsion  flow  in  a  porous 
medium  is  a  non-Darcy  flow  process,  because  the  degree  of 
pore  plugging  and  thus  the  permeability  of  the  medium  is 
dependent  on  the  pressure  gradient; 

2)  to  examine  the  flow  behaviour  of  an  emulsion 
through  various  porous  media; 
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3)  to  examine  the  change  in  the  droplet  size 
distribution  of  the  produced  emulsion  as  compared  to  the 
injected  emulsion  during  emulsion  flow; 

4)  to  develop  a  mathematical  model  that  describes 
emulsion  flow  in  a  porous  medium,  previously  saturated  with 
the  external  phase  of  the  emulsion.  This  model  assumes  that 
emulsion  droplets  partially  plug  off  the  porous  medium  and 
the  rate  of  plugging  can  be  modelled  by  a  rate  equation. 

The  coupling  of  this  kinetic  equation  of  plugging  with  a 
material  balance  equation  yields  a  model  which  predicts  the 
emulsion  oil  saturation  within  the  porous  medium.  This  type 
of  approach  has  been  successfully  used  for  the  flow  of 
suspensions  through  porous  media,  a  good  review  of  which  is 
given  by  Herzig,  Leclerc  and  LeGoff  (15)* 


CHAPTER  2 


THEORY  AND  LITERATURE  REVIEW 

2.1  EMULSION  RHEOLOGY 

In  this  study,  an  emulsion  viscosity 
correlation  was  used  to  evaluate  the  produced  emulsion  oil 
concentration, 

ve  =  i>w  (  i  +  c ,r  +  C2r2  )  (2.1) 

where 

ve  =  emulsion  viscosity  (  mm2  s"1  ) 

=  water  (  external  phase  )  viscosity  (  mm2  s"1  ) 

T  =  oil  (  internal  phase  )  concentration  (  -  ) 

Q ,  C2  =  constants  (  -  ) . 

The  factors  which  affect  emulsion  viscosity  are  reviewed 
briefly  and  the  incorporation  of  these  factors  in  the 
correlation  is  described.  The  review  was  based  on  the  books 
of  Becher  (16)  and  Sherman  (17)  and  therefore  these  books 
should  be  consulted  for  the  reference  to  the  original 
papers . 

Sherman  (18)  listed  six  factors  which  may  affect  the 
rheological  properties  of  an  emulsion: 

1)  Viscosity  of  the  external  phase  (  /xex  ), 

2)  Volume  concentration  of  the  disperse  phase  (  T  ), 

3)  Viscosity  of  the  internal  phase  (  fi\  )  t 

4)  Nature  of  the  emulsifying  agent  and  the 
interfacial  film  formed  at  the  interface, 

5)  Electroviscous  effects  and 
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6)  Particle  size  distribution. 

Each  of  these  variables  are  reviewed  individually  below. 

1)  Viscosity  of  the  External  Phase 

Virtually  all  emulsion  viscosity  equations  indicate  a 
direct  proportionality  between  the  external  phase  viscosity 
and  the  emulsion  viscosity.  The  general  form  of  the 
relationship  is: 

Me  =  Mex  .  f(x)  (2.2) 

where 

/ ie  =  emulsion  viscosity  (  Pa.s  ) 

Mex  =  viscosity  of  the  external  phase  (  Pa.s  ) 
f(x)  =  function  dependent  on  other  properties  which 
affect  the  viscosity  (  -  ) 

2)  Volume  Concentration  of  the  Internal  Phase 

The  principle  attention  of  researchers  in  the  field  of 
emulsion  viscosity  has  been  directed  towards  the  effect  of 
the  volume  concentration  of  the  internal  phase. 

The  classic  equation  relating  the  viscosity  of  a 
suspension  with  that  of  the  suspending  liquid  and  the 
volume  fraction  is  that  of  Einstein  (19) 

Me  =  /iw  (  l  +  2 . 5  r  )  (2.3) 

This  equation  is  probably  only  valid  for  values  of  r  less 
than  0.02.  Subsequent  workers  have  attempted  to  modify 
Equation  2.3  by  introducing  a  power  series  in  the  volume 
concentration.  These  take  the  general  form: 
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Me  -  Mw  (  i  +  C1  r  +  C2r  +  C3r  +  . . .  )  (2.4) 

where 

0^02,03=  constants 

and  C1  is  usually  taken  to  he  2.5,  as  in  Equation  2.3. 

Guth,  Gold  and  Simha  (20)  calculated  the  value  of  C2 
on  theoretical  grounds  and  arrived  at  the  equation 

Me  =  Mw  (  1  +  2 .5  r  +  14. ir2  )  (2.5) 

which  increases  the  validity  of  the  equation  to  a  value  of 
T  of  about  0.06. 

Oliver  and  Ward  (21),  analysing  data  on  model  emulsions 
consisting  of  suspensions  of  rigid  spheres  of  varying  size 
distributions,  discovered  that  a  straight  line  relationship 
existed  if  1  -  Mw/Me  was  plotted  versus  T  for  values  of 
the  volume  concentration  up  to  20  per  cent  (  or  slightly 
higher) .  Algebraically  expressed,  this  yields 

Me  =  Mw/  (  1  -  kr  )  (2.6) 

which  for  low  concentrations  can  be  written  as: 

Me  =  Mw  (  1  +  kr  +  k2r2  +  k3r3  + - )  (2.7) 

where  k  =  constant. 

Equation  2.7  is  a  form  of  Equation  2.4,  except  that  a 
functional  relationship  is  supposed  to  exist  between  the 
constants . 

For  emulsions  containing  more  than  5°  per  cent  of  the 
disperse  phase, it  is  generally  considered  they  will  show  a 
non-Newtonian  behaviour.  For  such  systems,  Hatschek  (22) 
derived  the  equation 

(<e  =  /*«/  (  1  -  r  3  ) 


(2.8) 
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Sibree  (23)  in  a  series  of  investigations  on  paraffin- 
in-water  emulsions,  found  it  necessary  to  modify  Hatschek's 
equation 

Me  =  Mw  /  |  1  -  (  hD  }  (2.9) 

where  the  constant  h  was  termed  the  volume  factor  by  Sibree . 
Sibree  found  that  h  was  nearly  1.3  for  a  number  of  oil-in¬ 
water  emulsions  of  a  variety  of  particle  sizes.  From  these 
data  it  was  concluded  that  the  droplets  behaved  as  though 
they  were  3°  percent  greater  in  volume  than  would  be 
expected  from  the  volume  concentration. 

The  majority  of  the  equations  found  in  the  literature, 
some  of  which  are  listed  above,  relating  emulsion  viscosity 
and  the  internal  phase  concentration,  under  certain 
conditions,  can  be  reduced  to  the  power  series  shown  as 
Equation  2.4. 

3)  Viscosity  of  the  Internal  Phase 

The  effect  of  the  viscosity  of  the  internal  phase  was 
introduced  to  correct  for  currents  set  up  within  the 
dispersed  droplets  under  the  shearing  flow  force,  under  the 
supposition  that  the  droplet  behaves  like  a  fluid  rather 
than  a  solid  particle .  This  effect  was  reported  by  Taylor 
(24),  who  extended  the  hydrodynamic  considerations  of  the 
classical  Einstein  relation,  Equation  2.2,  by  assuming  that 
any  interfacial  film  which  existed  merely  transmitted 
tangential  stress  from  one  phase  to  the  other,  and  arrived 
at  the  relation 
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fte  =  ft 


ex 


+  2.5  r 


Mi 


2/5  jJL 


ex 


Mi 


M, 


(2.10) 


ex 


which  can  be  rewritten  as 


Me  =  ft 


ex 


Mi 

i  +  r  5/2 


ex 


Mi_ 

Mex 


(2.11) 


where  /*.  =  viscosity  of  the  internal  phase . 

Therefore  according  to  Taylor's  equation,  the 
coefficient  in  Einstein's  equation  varies  from  1.0  to  2.5 
as  Mj /Mex  varies  from  zero  to  infinity.  The  effect  of 
increasing  the  internal  phase  viscosity  is  to  increase  the 
viscosity  of  the  emulsion. 

Taylor's  hydrodynamic  theory,  that  the  internal  phase 
viscosity  significantly  affects  the  emulsion  viscosity,  is 
valid  when  the  droplets  behave  as  liquids.  If  however,  the 
droplets  behave  as  rigid  spheres,  then  the  chemical  nature 
of  the  disperse  phase  has  a  greater  effect. 


4)  Interfacial  Film  and  Emulsifying  Agent 

The  interfacial  film  may  affect  the  emulsion  viscosity 
by  reason  of  its  effect  on  the  internal  circulation  and/or 
movement  of  the  droplet . 

Toms  (25)  made  an  intensive  investigation  on  a  series 
of  oil-in-water  emulsions  to  evaluate  the  effect  of 
changing  the  emulsifier  and  emulsifier  concentration  on  the 
emulsion  viscosity.  The  data  were  used  to  calculate  the 
volume  factor  h  of  Equation  2.9  and  the  average  value 
was  close  to  the  value  of  1*3*  Individual  variations  were 
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considerable  and  appeared  to  be  attributable  to  both  a 
variation  in  the  type  of  emulsifying  agent  and  the  internal 
phase .  Variations  in  the  emulsifier  concentration  appeared 
to  have  a  lesser  effect. 

5)  Electroviscous  Effect 

Smoluchowski  (26)  reported  that  particles  bearing  an 
electric  charge  showed  a  viscosity  exceeding  that  of  a 
similar  system  of  uncharged  particles.  The  contribution 
of  the  charge  was  to  increase  the  viscosity  of  the  emulsion, 
regardless  of  the  sign  of  the  charge .  Taking  this  effect 
into  account,  Smoluchowski  extended  Einstein's  equation  to 


where 

k  =  the  specific  conductivity  of  the  suspension 

_  i 

(  mhos  cm  ) 

€*  =  the  dielectric  constant  of  the  dispersion 

medium  (  coulomb  volt  cm  ) 

f  =  the  electrokinetic  potential  of  the  charged 

particles  (  volts  ) 

r,  =  the  radius  of  the  particles  (  cm  ) 
d 

6)  Particle  Size  Distribution 

Richardson  (27)  concluded  from  his  studies  that  the 
apparent  viscosity  of  emulsions,  having  the  same 
concentration  and  the  same  size  distribution  about  a  mean 
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droplet  diameter,  was  inversely  proportional  to  the  mean 
droplet  diameter.  The  observed  viscosity  variations,  of  a 
concentrated  emulsion,  with  concentration  and  rate  of  shear 
can  be  explained  in  terms  of  the  work  done  in  distorting  the 
droplets  and  sliding  them  past  each  other. 

Roscoe  (28)  obtained  two  modified  forms  of  the 
Einstein  equation,  one  of  which  applies  to  systems  in  which 
the  particle  size  distribution  is  broad 

Me  =  Mw  (  l  -  r  f5  (2.13) 

and  one  for  sharply-distributed  emulsions 

Me  =  M„  (  1  -  3.5  r  ]2'5  (2.14) 

Equation  2.13  is  applicable  over  all  concentration  ranges. 
Equation  2.14  is  appropriate  at  high  concentrations,  but  for 
concentrations  of  T  below  0.05  "the  Einstein  equation 
applies  to  such  sharply-distributed  systems. 

The  work  of  Sherman  (29)  showed,  that  for  water-in-oil 
emulsions  and  high  concentration  oil-in-water  emulsions, 
that  the  emulsion  viscosity  decreased  with  an  increase  in 
the  mean  droplet  size.  However,  for  lower  concentration 
oil-in-water  emulsions, the  emulsion  viscosity  was 
independent  of  the  droplet  size. 

In  summary,  the  factors  of: 

1)  the  viscosity  of  the  external  phase  and 

2)  the  volume  concentration  of  the  internal  phase  can 
be  accounted  for  with  an  equation  of  the  type 

Me  =  Mw  (  1  +  c,r  +  c/  ) 


(2.15) 
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The  factors  of 

3)  the  viscosity  of  the  internal  phase, 

4)  the  interfacial  film  and  emulsifying  agent, 

5)  the  electroviscous  effect  and 

6)  the  particle  size  distribution 

have  been  shown  to  affect  the  constants  in  the  Einstein  and 
power  series  models.  Therefore  on  defining  the  oil  and 
emulsifier  for  the  system,  the  effects  of  the  above  four 
factors  can  be  incorporated  into  the  constants  of  the  model. 
During  an  emulsion  flood  the  emulsion  parameters 

1)  emulsifier  concentration, 

2)  particle  size  distribution  and 

3)  emulsion  internal  phase  concentration 

will  vary  with  position  and  time.  The  first  two  parameters 
are  expected  to  have  only  a  small  effect  on  the  emulsion 
viscosity  for  a  low  concentration  (  T  =  0.15  )  oil-in¬ 
water  emulsion.  Therefore  a  correlation  of  the  type 

/*e  =  Mex  (  i  +  C,  r  +  C2r2  +  C3r3  )  (2.7) 

can  be  used  to  evaluate  the  produced  emulsion  concentration 
from  its  viscosity,  because  the  variable  which  has  the  most 
significant  effect  on  the  emulsion  viscosity  during  a  flood, 
is  the  concentration  of  the  disperse  phase . 

In  this  study,  the  correlation  was  defined  in  terms  of 
the  kinematic  rather  than  the  dynamic  viscosity  because  this 
was  the  parameter  measured.  This  also  removed  the  necessity 
of  measuring  the  density  of  the  produced  emulsion.  The 
general  model 
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"e  =  K  (  i  +  C,r  +  C2r2  +  c3r3  )  (2.16) 

was  compared  with  the  experimental  data  and  it  was  found 
that  Equation  2.1  could  provide  a  good  correlation. 

2.2  A  CONCEPT  OF  EMULSION  FLOW 

A  concept  required  in  dealing  with  emulsion  flow  in  a 
porous  medium  may  be  examined  by  considering  a  single 
droplet  of  oil  in  the  emulsion  entering  a  pore  constriction 
of  smaller  diameter,  schematically  illustrated  in  Figure 
2.1.  The  leading  edge  has  a  smaller  radius  than  that  portion 
of  the  droplet  still  in  the  pore,  thus  the  capillary 
pressure  is  greater  at  the  front  than  at  the  back  of  the 
droplet.  A  pressure  differential  is  required  to  force  the 
droplet  through  the  constriction.  The  pressure  differential 
is  given  by  the  equation 

AP  =  2  7  cos#  j  1  .  i 

\  r,  r> 

where 

AP  =  pressure  differential  (  Pa  ) 
y  =  interfacial  tension  between  the  oil  and  water 
(  N/m  ) 

0  =  contact  angle  (  -  ) 

r;  =  pore  throat  radius  (  m  ) 
r2  =  upstream  oil  droplet  radius  (  m  ) 

Thus  if  the  pressure  differential  across  the  droplet 
is  greater  than  AP,  then  the  oil  droplet  will  deform  and 
pass  through  the  constriction.  If  however,  the  pressure 
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FIGURE  2.1  SCHEMATIC  DIAGRAM  OF  AN 

EMULSION  OIL  DROPLET  TRAPPED 
BY  CAPILLARITY  AND  A  PORE 
CONSTRICTION 
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differential  is  less  than  AP,  then  the  oil  droplet  will 
plug  the  pore  constriction.  This  describes  the  basis  for 
permeability  reduction  of  the  porous  medium  by  emulsions. 

Conditions  conducive  to  pore  plugging  and  permeability 
reduction  can  be  implied  from  Equation  2.17  and  Darcy's 
Law 

q  =  _  K _ A  dP  (2.18) 

I u  dL 

where 

o 

q  =  volumetric  flowrate  (  cmJ  /  s  ) 

K  =  permeability  (  cm  ) 

A  =  cross-sectional  area  (  cm^  ) 

—  =  pressure  gradient  (  Pa  /  cm  ) 
dL 

/x  =  viscosity  (  Pa.s  ). 

These  conditions  are: 

1)  the  oil  droplets  within  the  emulsion  should  be 
larger  than  the  pore  throat  constrictions  in  the 
porous  medium, 

2)  a  higher  interfacial  tension  between  the  water  and 
oil  will  favour  the  blockage  of  the  pore 
constrictions , 

3)  the  flow  rate  has  to  be  less  than  some  critical 
value  in  order  that  the  pressure  differential 
acting  across  the  droplet  is  less  than  AP,  thus 
allowing  plugging  and 

4)  for  emulsions  with  the  same  maximum  droplet  sizes 
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and  interfacial  properties  and  flowing  at  the  same 
rate,  the  emulsion  with  the  lower  viscosity  will 
favour  plugging  because  the  pressure  gradient  will 
be  lower  for  such  systems. 

The  'non-Darcy'  flow  concept  of  emulsion  flow, 
described  by  McAuliffe  (12)  as  'pseudo  non-Newtonian',  may 
be  explained  by  examining  Equation  2.17.  At  lower  pressure 
gradients,  the  average  pore  size  which  can  be  plugged  will 
be  larger  and  the  pore  volume  plugged  or  isolated  within 
the  core  will  also  increase .  Due  to  the  increase  in  the 
pore  volume  unavailable  for  fluid  flow,  the  permeability 
of  the  core  will  be  reduced;  thus  the  core  permeability 
is  a  function  of  the  pressure  gradient.  The  dependence 
of  the  permeability  on  the  pressure  gradient  makes  emulsion 
flow  a  'non-Darcy'  flow  process. 

A  topic  allied  to  emulsion  flow  is  the  mobilisation 
of  ganglia  and  its  influence  on  the  oil  recovery  process. 
Both  processes  involve  understanding  oil  droplets  and 
their  flow  behaviour  through  constrictions.  Two  of  the 
critical  properties  are  the  maximum  oil  droplet  size  and 
the  minimum  pore  constriction  size.  These  dimensions 
determine  the  maximum  pressure  gradient  necessary  to 
mobilise  an  oil  ganglia  and  the  maximum  pressure  gradient 
below  which  plugging  will  occur  in  emulsion  flow.  More 
work  has  been  carried  out  on  the  mobilisation  of  ganglia 
(3°)  ~  (35)  than  on  emulsion  flow.  A  useful  concept 
relating  the  capillary  number  necessary  to  mobilise  oil 
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droplets,  has  been  developed  (3°)  -  (3*0*  The  application 
of  this  type  of  concept  for  describing  emulsion  flow  is 
discussed  in  Section  5*3  of  this  study. 


2.3  REVIEW  OF  EMULSION  FLOW  CHARACTERISTICS 

McAuliffe's  laboratory  work  (12)  showed  that  oil-in- 
water  emulsions  could  effectively  reduce  the  water 
permeability  of  sandstone  cores  when  the  initial  water 
permeability  was  less  than  two  darcies.  The  permeability 
reduction  was  retained  even  when  the  flowing  emulsion  was 
followed  by  many  pore  volumes  of  water. 

McAuliffe  further  showed  that  the  degree  of 
permeability  reduction  increased  with  an  increase  in  oil 
droplet  size  and  decreased  with  an  increase  in  the  pressure 
differential.  He  described  the  flow  of  oil-in-water 
emulsions  through  a  porous  medium  as  'pseudo  non-Newtonian' , 
regardless  of  the  emulsion  concentration. 

Results  on  the  flow  of  an  oil-in-water  emulsion 
through  three  cores  of  different  permeabilities  (  1460,  595 
and  395  md  )  mounted  in  parallel,  showed  that  the  core 
permeabilities  were  reduced  by  factors  of  35 »  21  and  12 
times,  respectively.  The  results  of  this  experiment  also 
showed  that  initially  a  large  amount  of  emulsion  entered 
the  more  permeable  core.  As  this  occurred,  the  flow  became 
more  restricted  and  the  emulsion  began  to  flow  into  the 
less  permeable  cores.  Eventually  the  same  number  of  pore 
volumes  passed  through  each  core,  resulting  in  an 


. 
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improvement  of  the  displacement  or  pore  volume  sweep 
efficiency.  A  further  observation  of  this  experiment  was 
that,  apparently,  emulsified  oil  is  much  more  efficient  at 
reducing  the  permeability  of  the  core  than  is  oil  left 
after  waterflooding. 

McAuliffe 's  results  also  showed  that  oil-in-water 
emulsions  displaced  oil  from  sandstone  core  material  more 
efficiently  than  did  water  alone.  McAuliffe  believed  this 
was  due  to  a  correction  of  heterogeneity  of  the  porous 
medium . 

McAuliffe  (10)  reported  the  results  of  an  emulsion 
flood  field  test  which  was  conducted  to  see  if  oil-in¬ 
water  emulsions  would  reduce  channeling  from  injection  to 
producing  wells.  A  conclusion  of  the  field  test  was  that 
the  emulsified  oil  decreased  fingering  and  increased 
volumetric  sweep  efficiency.  These  results  were  based  on 
the  following  observations: 

1)  increased  oil  production  and  lower  WOR's  from  the 
wells  surrounding  the  emulsion  treated  injectors, 
compared  with  little  or  no  increase  in  oil 
production  and  increasing  WOR's  in  the 
surrounding  conventional  waterflooding  patterns, 

2)  production  of  a  greater  proportion  of  saline 
formation  water  from  several  wells  after  emulsion 
injection,  indicating  a  decrease  in  the  injection- 
water  channeling  and 

3)  a  favourable  change  in  the  flood  pattern  after 
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emulsion  injection,  as  suggested  by  radioactive 
tracer  analysis. 

The  work  of  Uzoigwe  and  Marsden  (36)  showed  that  no 
flow  restriction  occurred  when  an  oil-in-water  emulsion  was 
flowed  through  glass  bead  packs.  This  was  probably  because 
the  pore  throats  in  the  unconsolidated  glass  bead  pack 
were  larger  than  the  maximum  size  of  the  oil  droplets. 

Alvarado  (37)  reported  that  for  the  conditions  of  his 
study,  oil-in-water  emulsions  would  reduce  the  core 
permeability.  Over  the  flow  rates  investigated,  312  -  1116 
ft/day,  the  reduction  was  independent  of  flow  rate  and 
pressure  gradients  and  thus  the  flow  of  emulsions  was 
classified  as  a  Darcy  flow  process.  At  these  rates,  the 
pressure  drops  varied  from  2 6  -  150  psi/ft.  Both  the 
flow  rates  and  the  pressure  drops  were  much  higher  than 
those  at  normal  field  conditions,  which  are  of  the  order  of 
1  ft/day  and  1-5  psi/ft. 

Alvarado  further  reported  that  permeability  reduction 
increased  with  emulsion  concentration,  which  conflicts  with 
the  results  described  in  Section  2.2.  Permeability  reduction 
should  increase  with  a  decrease  in  the  emulsion  viscosity 
and  thus  a  decrease  in  the  concentration.  However, 

Alvarado's  observation  was  based  on  studies  where  the  core 
permeability  was  different  for  each  experiment.  The  core 
permeability  was  higher  for  higher  emulsion  concentrations. 

As  McAuliffe  has  already  shown  that  permeability  reduction 
is  a  function  of  the  core  permeability,  it  is  believed  that 


Alvarado's  observation  is  invalid  as  two  parameters  were 
varied  simultaneously. 

Cartmill  (38)  carried  out  experiments  on  the  flow  of 
oil-in-water  emulsions  to  establish  if  oil  migrates 
through  the  reservoir  sands  in  the  form  of  an  emulsion.  He 
found  retention  of  the  oil  at  the  interface  between  the 
coarse  and  fine  media  in  the  core  packs,  and  as  a 
consequence,  the  permeability  of  the  porous  medium  was 
reduced.  Since  the  size  difference  between  the  lj-  micron 
oil  droplets  and  the  porous  media  of  37  -  88  microns  was 
substantial,  he  reasoned  that  electrostatic,  rather  than 
capillary  forces,  were  the  cause  of  the  oil  accumulation. 

Devereux  (39)  showed  that  the  amount  of  oil  retained 
in  a  porous  medium  during  emulsion  flooding  decreased  as  th 
pressure  differential  increased.  Also  over  the  range  of 
emulsion  concentrations  (  0.001  -  0.05  )  and  pressure 
differentials  (  8  -  100  psi/ft  )  studied,  the  retained  oil 
concentration  increased  with  an  increase  in  the  inlet 
emulsion  concentration  for  the  pressure  range  from  8  -  40 
psi/ft.  For  the  pressure  range  from  60  -  100  psi/ft,  the 
retained  oil  concentration  passed  through  a  maximum  at  an 
emulsion  concentration  of  0.01. 

The  results  of  Devereux  also  showed  that  for,  low  oil 
concentration  emulsions,  the  oil  saturation  front  moved 
through  the  core  at  a  slower  velocity  than  the  injection 
velocity,  because  the  oil  saturation  within  the  invaded 
region  was  much  higher  than  the  injected  emulsion 
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concentration . 

The  general  flow  characteristics  of  emulsion  flow  are 
reviewed  above,  the  consequences  of  which  are  listed  below: 

1)  emulsions  are  more  efficient  at  displacing  oil 
than  either  water,  caustic  or  acid  displacing 
agents , 

2)  emulsion  flooding  improves  the  volumetric  sweep 
efficiency  and  decreases  fingering  (10)  and 

3)  the  'non-Darcy*  flow  nature  is  an  advantage  when 
an  emulsion  is  injected  into  an  oil  reservoir, 
where  the  largest  change  in  pressure  for  a  unit 
distance  occurs  near  the  wellbore.  As  the  pressure 
drop  per  unit  distance  decreases  logarithmically 
with  distance  from  the  injection  well,  emulsions 
can  be  injected  into  the  formation  easily  and  will 
be  most  effective  in  reducing  the  permeability  at 
some  distance  away  from  the  wellbore  (10). 


2.4  EMULSION  FLOW  MODELS 

Two  models  for  emulsion  flow  in  a  porous  medium  were 
found  in  the  literature.  The  model  of  Spielman  and  Goren 
(40)  is  related  to  the  coalescence  and  flow  of  an  emulsion 
in  a  porous  bed,  and  is  specifically  related  to  oil-water 
separators.  As  such  it  will  not  be  discussed  further  in  this 
thesis . 

The  model  of  Devereux  (  4l )  is  related  to  the  flow  of 
an  oil-in-water  emulsion  in  a  porous  medium  previously 
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saturated  with  water.  Devereux  derived  an  equation  for  the 
one-dimensional  flow  of  an  emulsion  in  a  porous  medium  and 
the  relationship  found  is  as  follows: 


t 

L  </>  S 


M  -  Mi  -  Me)  Q  T 
Mj  (M  -  Me  )7J  L  </>  SQ 
M  M-,  1'  +  (  M  -  M;  -  Me  )  MP/L 


In 


Mj  (  M  -  Me  )27'2  j 

-  (  M  -  Me  )T'  Q  r  } 
MP/L  -  M  y  L  <t>  S0  ) 


(2.19) 


where 

t  =  time  (  seconds  ) 

L  =  length  of  solid  (  cm  ) 

9  =  porosity  (  -  ) 

SQ  =  oil  saturation  in  the  invaded  region  (  -  ) 

M  =  absolute  mobility  (  permeability  /  viscosity  ) 
of  the  porous  media  to  the  external  phase  with 
no  internal  phase  in  the  porous  media  (  cm2/ 
poise  ) 

Mj  =  mobility  of  the  internal  phase  during  the 
emulsion  flood  (  cm2/  poise  ) 

Me=  mobility  of  the  external  phase  during  the 
emulsion  flood  (  cm2/  poise  ) 

/ 

^  =  capillary  constant  (  dyn  /  cm  ) 
r  =  emulsion  concentration,  volume  of  internal 
phase  to  volume  of  emulsion  (  -  ) 

Q  =  cumulative  injected  emulsion  per  unit  cross- 
sectional  area  (  cm3/  cm2  ) 
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P  =  pressure  differential  (dyn  /  cm  ) 

The  pertinent  parameters  of  this  model  are  shown  in 
Figure  2.2.  Within  the  invaded  region,  the  internal  phase 
flows  at  a  rate  q; ,  with  mobility  Mj  and  the  external  phase 
flows  at  a  rate  qe,  with  mobility  Me  .  Beyond  this  region, 
the  external  phase  flows  at  a  rate  q.  The  saturation  is 
presumed  to  have  a  uniform  value  of  SD  in  the  invaded 
region,  changing  discontinuously  to  zero  in  the  uninvaded 
region.  The  model  assumes  incompressibility  and  absence  of 
gravitational  effects  but  capillary  effects  were  not 
neglected . 

The  driving  force  for  flow  of  the  internal  phase  was 
assumed  to  consist  of  the  pressure  gradient  in  the  external 
phase  minus  a  retarding  force  due  to  capillary  effects.  The 
retarding  force  was  attributed  to  three  phase  (  internal 
phase,  external  phase  and  solid  )  contact  and  this  may 
represent : 

a)  contact  angle  hysteresis  for  an  emulsion  droplet 
moving  through  an  essentially  cylindrical 
capillary, 

b)  variation  in  the  length  of  advancing  versus 
receding  contact  line  for  a  droplet  moving  through 
a  pore  of  varying  cross  section  with  whose  walls  it 
forms  a  well-defined  contact  angle, or 

c)  increase  in  surface  area  due  to  deformation  of  a 
drop  moving  through  a  pore  of  varying  cross 
section  with  which  it  forms  no  contact. 


. 
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L  ,  P 


V 


FIGURE  2.2  DEFINITION  OF  THE  FLOW 

PARAMETERS  FOR  LINEAR 
FLOW  GEOMETRY  MODEL 
OF  DEVEREUX 
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These  three  effects  are  shown  schematically  in  Figure  2.3* 
The  capillary  factor  is  thus  a  summation  of  several  effects 
which  cannot  he  expressly  calculated.  Simple  expressions 
can  he  used  to  calculatedthe  typical  contrihutions  of  the 
suggested  models  to  the  capillary  factor.  The  capillary 
force  due  to  contact  angle  hysteresis  is  given  hy 

l'  =  ^  ^o  d  ^  (  cos  ^  -  cos  02  )  (2.20) 

dd 

where 

y  =  surface  tension  (  dyn  /  cm  ) 
d'  =  pore  diameter  (  cm  ) 
d^  =  droplet  diameter  (  cm  ) 

6  =  contact  angle  (  -  ) 

1  =  refers  to  leading  edge  of  the  droplet 

2  =  refers  to  trailing  edge  of  the  droplet 

The  assumptions  made  in  the  derivation  of  the  model  were: 

1)  the  emulsion  does  not  coalesce,  or  adhere  to  the 
solid  surface, 

2)  ]VK,  and  7/  are  functions  only  of  saturation  SQ , 

3)  the  saturation  is  invarient  in  time  and 

4)  the  pressure  differential  across  the  system  was 
invarient  in  time . 

The  consequences  of  the  assumptions  described  above  are: 

1)  that  plugging  of  the  medium  is  not  incorporated  in 
the  model,  as  all  of  the  oil  phase  is  assumed  to  he 
mobile , 
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b)  Contact  Line  Variation 


FIGURE  2.3  POSSIBLE  MECHANISMS 

LEADING  TO  A  CAPILLARY 
FACTOR 


27 


2)  the  oil  saturation  front  is  treated  as  a  moving 
boundary,  rather  than  as  a  front  with  a  saturation 
distribution . 

With  all  the  physical  data  known  or  recorded,  the  three 
constants  M^,  and  l'  were  then  fitted  to  the 
experimental  data  using  a  least  squares  method.  The  flow 
model  was  a  good  exhibit  of  the  results  for  high  pressure, 
low  emulsion  concentration  systems.  The  model  did  not  fit 
the  behaviour  for  a  high  emulsion  concentration,  low 
pressure  system.  An  explanation  for  the  poor  description 
of  the  high  concentration,  low  pressure  systems  was  the 
violation  of  assumptions  1  and  3  above.  Both  McAuliffe  (12) 
and  Jennings  et  al  (3)  have  reported  that  droplets  become 
lodged  in  pore  constrictions,  which  conflicts  with 
assumption  1.  While  Devereux  reported  that  the  saturation 
is  time  dependent  and  therefore  conflicts  with  assumption 

3- 

Further,  Devereux  assumed  that  all  the  oil  in  the 
invaded  region  is  mobile  and  contributes  to  the  flow 
constant.  It  is  more  likely  that  the  majority  of  this  oil 
is  lodged  within  a  pore  constriction  and  only  a  small 
proportion  of  the  oil  is  flowing  within  the  core. 

Thus  the  conditions  under  which  Devereux' s  model  provide 
provides  a  good  correlation,  are  the  conditions  most 
unlikely  to  be  found  in  an  actual  reservoir.  Therefore, 
another  concept  for  emulsion  flow  needs  to  be  formulated. 


CHAPTER  3 

EMULSION  FLOW  MODEL 

3.1  INTRODUCTION 

At  present,  immiscible  displacement  theory 
is  based  upon  Darcy's  Law  for  multiphase  flow,  an  equation 
defining  capillary  pressure  and  the  conservation  of  mass  in 
the  porous  medium.  This  formulation  was  originally  derived 
for  the  non-capillary,  two  phase  flow  model  of  waterflooding 
by  Buckley  and  Leverett  (42)  and  has  since  been  extended  to 
predict  other  flooding  processes  such  as  miscible  flooding 

(43). 

Inherent  in  the  Buckley-Leverett  solution  and 
subsequent  numerical  solutions  of  the  immiscible 
displacement  equation,  is  that  one  or  both  of  the  effective 
permeability  -  saturation  and  the  capillary  pressure  - 
saturation  relationships  must  be  known.  This  formulation 
requires  experimentally  determined  properties  and/or 
mathematical  models  for  the  above  relationships.  Due  to 
these  inherent  uncertainties  and  the  problem  of  defining  the 
fractional  flow  equation  for  emulsion  flow,  it  was  felt  that 
a  further  extension  of  their  theory  to  predict  emulsion  flow 
was  not  feasible  and  a  new  concept  for  emulsion  flow  would 
need  to  be  formulated. 

For  this  new  theory  it  was  assumed  that  the  plugging 
of  the  porous  medium  by  oil  droplets  could  occur  and  that 
the  rate  of  plugging  could  be  modelled  by  a  rate  equation. 
The  combination  of  this  rate  equation  and  a  material  balance 
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equation,  was  analytically  solved  and  resulted  in  a  one- 
dimensional  saturation  profile  for  emulsion  flow  in  a 
porous  medium,  previously  saturated  with  the  external  phase 
of  the  emulsion. 

3.2  MATERIAL  BALANCE  EQUATION 

It  is  assumed  that  both  the  water  and  oil  phases  are 
incompressible  and  therefore  a  volume  balance,  rather  than 
a  mass  balance,  could  be  carried  out  over  the  core. 

Having  assumed  that  oil  droplets  can  plug  off  pore 
constrictions,  then  two  types  of  oil  might  be  assumed  to 
exist  in  the  core : 

1)  stationary  oil,  plugged  within  pore  constrictions 
and 

2)  moving  oil,  within  the  flowing  emulsion. 

Considering  a  porous  medium  of  length  AL,  area  A  and 
porosity  0,  the  expression 

A  AL0a  represents  the  volume  of  retained  or  plugged 
oil  and 

A  AL0er  represents  the  volume  of  flowing  oil  within 
the  emulsion, 

where 

o  =  volume  of  retained  (plugged)  oil  per  unit  pore 
volume  (  -  ) 

T  =  emulsion  concentration  (  -  ) 

0e  =  porosity  through  which  the  emulsion  is  flowing 
(  -  ) 
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Associated  with  the  retained,  immobile  oil  droplets,  is  a 
volume  of  immobile  water  given  by 
0  A  AL  (  b  -  1  )  a 


where  b  =  a  blockage  coefficient  defined  as 

volume  immobile  oil  +  volume  immobile  water 
volume  of  immobile  oil 


such  that 


0  =  </>e  +  b  O 

The  general  volume  material  balance  is  given  by 


A  (  0a  +  0er  ) 

+ 

a 

qr 

ai_ 

(3-D 


(3-2) 


where  q  =  volumetric  emulsion  flow  rate  (  cc/hr  ) . 

For  a  constant  flow  rate  experiment,  Equation  3*2 
simplifies  to 

d[a  +  t  r  ]  +  1  dr  =  0  (3-3) 

3t  a  <t>  JT 


which  can  be  written 

as0  +  q  ar  =0  (3-4) 

at  a  0  aL 

where  e  =  0e/0 

SQ =  oil  saturation,  total  volume  of  oil  per  unit 
pore  volume  (  -  ) 


3.3  PLUGGING  MECHANISMS  -  THE  KINETIC  EQUATION 

Unlike  the  material  balance  equation,  the  form  of  the 
rate  equation  is  unknown  and  has  to  be  derived  from  actual 
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experimental  data.  Before  discussing  the  possible  forms  of 
the  rate  equation,  the  various  oil  retention  sites  and 
capture  mechanisms  will  be  considered. 

The  most  likely  retention  site  for  larger  particles, 
described  in  Section  2.2,  is  the  constriction  site,  where 
the  droplet  is  captured  by  capillarity  when  the  pressure 
differential  is  not  sufficient  to  squeeze  the  droplet 
through  the  pore  constriction.  Another  likely  retention 
site,  is  a  dead  end  cavern  site  behind  a  plugged 
constriction  site . 

The  results  of  Cartmill  (38)  showed  that  small 
particles  tend  to  collect  on  the  solid  surface  and  to  each 
other,  whenever  the  pore  diameter  of  the  porous  medium  was 
abruptly  reduced.  The  suggested  reason  for  the  oil 
accumulation  was  electrostatic  forces  rather  than  capillary 
forces.  These  long  chains  and  clusters  of  droplets,  formed 
by  agglomeration  and  flocculation,  were  found  on  the 
surface  and  crevice  sites.  A  schematic  diagram  showing  each 
of  these  retention  sites  is  given  in  Figure  3*1* 

In  emulsion  flooding,  as  in  waterflooding,  the  flood 
can  be  divided  into  three  zones: 

zone  1  :  behind  the  front,  where  the  retained  oil 
saturation  has  reached  its  maximum  value  crm,  and  the 
emulsion  is  flowing  with  the  inlet  concentration  F,  .  In 
this  zone  the  pressure  gradient  is  constant  in  a  linear 
flow  system  and  the  highest  in  the  emulsion  flow  system. 
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Surface  Site  Crevice  Site 


FIGURE  3.1  RETENTION  SITES 
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zone  2  :  in  the  front,  where  the  oil  retained  is 
between  zero  and  its  maximum  saturation  and  the  flowing 
emulsion  concentration  is  between  zero  and  the  inlet 
concentration  Tj  .  In  this  zone  the  pressure  gradient 
increases  with  distance  behind  the  leading  edge  of  the 
front . 

zone  3  s  ahead  of  the  front,  where  only  water  is 
flowing.  In  this  zone  the  pressure  gradient  is  constant  in 
a  linear  flow  system  and  the  lowest  in  the  emulsion  flow 
system . 

Considering  a  point  in  time,  a  droplet  may  become 
constricted  in  a  pore  throat,  then  as  zone  2  continues  to 
pass  through  this  position,  the  pressure  gradient  may 
increase  sufficiently  to  deform  the  oil  droplet  and  cause 
the  breaking  away  of  the  droplet.  Therefore  within  zone  2, 
there  will  be  two  processes: 

1)  plugging  of  the  medium  by  the  oil  droplets  and 

2)  the  possible  breaking  away  of  the  oil  droplets 
from  their  capture  sites. 

The  droplets  that  break  away  then  flow  through  the  medium, 
where  they  may  find  another  pore  constriction  small  enough 
for  capture .  Thus  within  the  system  there  are  two  types  of 
capture : 

l)  permanent  capture  :  where  the  droplet  becomes 
constricted  and  stays  there  while  the  system 
flow  rate  remains  unchanged  and 
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2)  semi-permanent  capture  :  where  while  the  pressure 
gradient  is  less  than  the  critical,  the  droplet 
remains  constricted.  However  when  the  pressure 
gradient  reaches  a  critical  level,  the  oil  droplet 
breaks  away  from  the  capture  site . 

Gruesbeck  and  Collins'  data  (44)  on  the  entrainment  and 
deposition  of  fine  materials  in  a  porous  media 
substantiated  that  the  above  phenomenon  does  occur. 

Modelling  the  system  of  permanent  and  semi-permanent 
capture  would  be  extremely  difficult.  Therefore  the  rate 
equation  proposed  in  the  mathematical  model,  describes  the 
rate  at  which  droplets  find  permanent  capture  sites,  with 
the  breaking  away  of  the  droplets  not  being  considered. 

This  is  the  approach  taken  by  Gruesbeck  and  Collins, 
although  they  reported  that  they  could  modify  their  theory 
to  include  this  stochastic  effect. 


Although  a  number  of  different  retention  sites  and 
mechanisms  for  capture  may  occur  within  the  medium,  it  is 
assumed  that  the  fractional  retention  can  be  defined  by 
the  equation 

f  =  A  </>  Al  da  1  (3-5) 

dt  q  r 

and  the  fractional  retention  per  unit  length  as 


A  0  da 
q  r  dt 


(3.6) 


where  the  fractional  retention  is  defined  as  the  volume  of 


oil  which  becomes  plugged  in  pore  constrictions  as  a 
fraction  of  the  total  oil  which  entered  an  element  of  the 
porous  medium. 
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The  above  equation  can  be  defined  in  terms  of  unit 
time  rather  than  length,  by  introducing  the  time  interval 
for  passage  through  the  unit  length,  which  is  q  /  A </)e. 

Thus 

frt  =  frl  <3-7) 

A  4>e 

which  after  substituting  into  Equation  3.6  yields: 

f„i  =  — —  0.8) 

rl  er 

Equations  3*6  and  3-8  can  be  rearranged  to  yield 


do 

at 


(3.8) 


do  q  r 

at  rl  X0 


(3.6) 


where  f  ^  and  f  ^  depend  on  the  parameters  which 
characterise  the  emulsion  flow  system  and  particularly  on 
the  retention  parameter  o  .  Differential  Equations  3*8  and 
3.8  are  analogous  to  the  classical  equations  of  chemical 
kinetics  and  for  this  reason  they  can  be  called  the 
kinetic  equations  of  plugging.  In  fact,  they  form  one 
kinetic  equation  only: 


at 


(3.9) 


where  H  is  called  the  plugging  rate  constant  and  has  the 


36 


-1 

dimensions  of  time  and  is  a  function  of  the  retained 
oil  saturation.  The  plugging  rate  constant  is  related  to 
f  ^  and  f  ^  by  the  equation 


H  = 


q. 

A  0 


(3-10) 


The  parameters  which  have  an  effect  on  the  clogging 
kinetics  are: 

1)  the  emulsion  which  is  characterised  by  its 
viscosity  fie ,  flow  rate  q,  the  droplet  size  d^, 
and  the  concentration  r,  and 

2)  the  porous  medium  which  is  characterised  by  its 
permeability  K,  porosity  0 ,  and  the  grain  size 

V  As  the  bed  clogs>  the  parameters  K-  and 

H  vary  with  the  oil  retention  a. 

Therefore  in  general,  the  rate  constant  H  can  be  written 
H  =  function  ( fie,  o  ,  0  ,  K  ,  q  ,  d  ,  d^  ,  0e  ) 
and  more  specifically 

H  =  HQ  f(  <7  )  (3.11) 

where 

Hq  =  function  (  /*e>  0  ,  K  ,  q  ,  d  ,  dd  ) 
and  the  parameters  of  a  and  0e  are  incorporated  into  the 
function  f  (  cr  )  . 

The  form  of  f(  a  )  is  unknown  and  must  be  assumed  and 
the  resultant  emulsion  flow  model  compared  with  the 
experimental  data.  The  three  forms  of  f(  a  )  investigated 
in  this  study  are 

case  1:  f(a)  =  (om-<r)  (3*12) 


case  2  : 

f  (  *  )  = 

(  °m  " 

O'  ) 

( 

/3  -  o-  ) 

(3-13) 

case  3  : 

f(  o  )  = 

(  °m  ~ 

O'  ) 

( 

/3  +  o-  ) 

(3-14) 

where  0  is  a  constant  for  the  kinetic  equation. 

An  illustration  of  the  forms  of  the  three  rate 
Equations  3*12  -  3*14  are  shown  in  Figure  3-2.  For  the 
kinetic  Equations  3*12,  3*13  and  3-14  when  am<  (3  ,  the 
plugging  rate  is  highest  at  time  t  =  0,  when,  as  yet,  no 
plugging  has  occurred.  With  the  growth  of  the  plugged 
constrictions,  the  rate  decreases,  reaching  a  zero  value 
after  some  time  t;  at  this  point  plugging  is  complete. 

The  kinetic  Equation  3*14,  when  cm>  (3  ,  assumes  that 
the  surface  activity  of  the  medium  for  droplet  capture  is 
greater  when  some  plugging  has  occurred.  An  explanation 
for  this  form  of  kinetic  equation  is  that  cavern  sites 
and  dead  zones,  as  well  as  constriction  sites,  have  been 
created,  thus  increasing  the  number  of  capture  sites. 
However,  the  increase  is  not  permanent,  with  eventually  the 
rate  decreasing  due  to  the  number  of  sites  available,  with 
the  rate  reaching  zero  after  some  time  t. 

3.4  INTEGRATION  OF  THE  DIFFERENTIAL  EQUATIONS 

Although  the  three  kinetic  Equations  3*12  -  3*14  are 
compared  in  this  study,  the  integration  of  the  differential 
equations  will  only  be  illustrated  using  Equation  3-14.  Only 
the  solutions  and  assumptions  made  in  solving  Equations 
3.12  and  3*13  are  reported. 

The  equations  to  be  solved  for  the  unknowns  r(L,t) 


. 


Plugging  Rate  do/d t 


i 


FIGURE  3.2  ILLUSTRATION  OF  THE  FORMS 

OF  THE  THREE  KINETIC 
EQUATIONS 
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and  a(L,t)  are 
do-(L,  t ) 


=  H  (  <Jm  -  o-(L,t))  (  0  +  a(L,t))  r(L,t) 
dt  u 


(3-15) 


and 

d[a(L,t)  +  er(L,t)] 
dt 


A  </> 


dr(L,t) 

dL 


=  o 


(3.2) 


3.4.1  Plugging  Process  at  the  Inlet  End  of  the  Porous  Medium 
The  form  of  the  plugged  oil  saturation  at  the  core 
inlet  end  (  L  =  0  )  can  be  determined  from  Equation  3.14. 

For  a  constant  flow  rate  case,  with  constant  emulsion 
injection  concentration  T- 

r(0, t )  =  r  =  constant  (3*16) 

and  for  the  porous  medium  initially  saturated  with  water 

*(0,0)  =  0  (3-17) 
then  Equation  3*14  can  be  written  as 
0, t ) 


dt 


=  Hq  (  <rm  -  cr(  0 ,  t )  )  (  0  +  <r(0,t))r. 


from  which  the  inlet  end  plugged  oil  saturation  with  time 
is  established: 


*  ( 0 ,  t )  = 


exp  (  H0t  P  (  <xm  +  0  ) )  -  1 


'm 


(3-18) 


exp  (  H0t  T;  (  +  &  )  )  + 

A  graphical  illustration  of  Equation  3*18  is  shown  in  Figure 
3-3*  The  function  demonstrates  that  the  oil  saturation  at 
the  inlet  end  of  the  core  does  not  jump  instantaneously  to 
its  maximum  value  when  flow  begins,  but  takes  some  time  to 


reach  this  value . 


. 


FIGURE  3.3  THE  PLUGGING  PROCESS  AT  THE 

INLET  END  OF  THE  POROUS 
MEDIUM 
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3*4.2  Solution  of  the  Retained  Oil  Saturation  Equation 


Eliminating,  in  Equation  3.14,  the  function  r(L,t) 
and  placing  the  parameter  <r(L,t)  in  a  differential 
operator  one  obtains 


r(L,t)  =  -  hc 


1 


—  Om 

0-+  /3 


(3.19) 


and  differentiating  Equation  3.19  with  respect  to  time 
produces 


dr(L,t) 

at 


[  v-  0] 


o'  +  j8 


(3-20) 


and  with  respect  to  length 


dr(L,t)  _  1 

dL  0  P  +0] 

Simplification  1: 

Now  in  order  to  solve  Equations  3.3  and  3-14, 
Equation  3*3  has  to  be  approximated  to 


at  ai_ 


In 


a  —  a. 


JD 


o'  +  (3 


(3-21) 


A  _i_  7  ar  t  g  ar  ^ 

at  at  a  0  aL 


(3.22) 


where  e  is  the  average  value  for  the  fraction  of  the  pore 
space  occupied  by  the  flowing  emulsion  within  the  emulsion 
front.  Without  this  simplification  the  equation  cannot  be 
solved  analytically  because  in  the  function  a  +  eF ,  e  is 
a  function  of  a 

e  =  1  -  b  a  (3-23) 

The  relationship  between  6  ,  b  and  a  is  unknown.  Thus  the 
function  a  +  eT  is  replaced  by  a  + t F  and  a  solution  of 


■ 


* 
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the  equation  is  possible . 


Introducing  the  definition 


O  -  <J 


m 


a  +  0 


=  G 


(3-24) 


and  substituting  Equations  3.24,  3.20  and  3.21  into  3.22 
one  obtains 


—  ^-2-  In  G  +  e  InG  ~  H0(am+ 0)  d 


A  </>  dt  3L 


dt 


at 


Oin+  /3G 


1  -G 


-  o  (3-25) 


Integration  of  Equation  3*25  with  respect  to  time  yields 


q  JLlnG  +  e  i_  InG  -H0(*m+/3) 

a  <i>  <3l  at 


°’m  +  ffG 

1  -  G 


=  C4  (3-26) 


The  constant  of  integration  is  evaluated  using  the 
following  boundary  condition 


<UL,t)  =0  t  <  — —  (3.27) 

vf 

where  v^  =  q  /  A  </>g  =  frontal  velocity. 

This  condition  states  that  there  is  no  oil  in  the  core  at 
the  position  L  before  the  emulsion  front  has  arrived. 

Thus  for  the  boundary  condition  3*27,  equals  zero 
and  Equation  3*26  becomes 

Q  a  InG  +  6  a  InG  ~  HA(  °m  +  0  )  f^m  +  /3Gl  _  n  /  o  oq  n 

a 0  aL  aT  L  i -g  J  13  J 

A  solution  to  Equation  3.28  in  this  form  is  possible, 
however  the  resulting  equation  is  too  complicated  for 
practical  purposes.  Thus  a  simplified  solution  will  be 
presented  which  yields  a  more  manageable  expression. 


' 
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Simplification  2: 

By  linearising  the  third  term  in  Equation  3*28 
around  <r(L,t  )  =  0  and  by  introducing  the  new  variable 


t  —  t  -  L 


Vr 


then  Equation  3-28  becomes 

_£L_  dG(  L  t)  -  H  p2 ( G  |  ffm  \ G 
A0  5L  \  4 

The  integration  of  Equation  3*3°  yields 


=  0 


G(L,t)  =  - 


m 


/3  ~  C5  exp  {- 


-,-1 


q 


(3-29) 


(3-30) 


(3.31) 


The  constant  of  integration  C ^  can  be  found  from  Equation 
3*18.  At  L=  0,  Equation  3*31  becomes 


(  L ,  r)  =  <%  (  C,  -  0  ) 


-1 


(302) 


where  from 


c5  =  0 


1  -  exp]  Hor\  (%+  0)t 


(303) 


On  substitution  of  the  expression  for  into  Equation  3-31 
one  obtains 


cr  (  L  ,  T  )  —  °m 


exp  td  -  1 


exp  + 


1  +  q~m 

4 


exp  1^  - 1 


(304) 


ld  =  H0I|  (  "m  +  0  )  T 

1  _  H0L A<£ 

d  q 


where 
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3*4.3  Solution  in  Terms  of  the  Flowing  Emulsion  Concentration 


Substituting  for  a(L,t)  and  its  derivative  from 
Equation  3*3^  into  Equation  3*15  and  solving  yields 


r(L,r)  =  r.  - — p-td - 

1  exp  td  +  exp  1^-1 


(3-35) 


Thus  Equations  3*24  and  3*25  are  the  solutions  of  the  system 
of  Equations  3*3  and  3*15  using  the  boundary  conditions  3*16, 
3.17,  3*2?  and  3*18  with  simplifications  1  and  2. 


3.4.4  Comparisons  of  the  Flow  Models 


The  solutions  using  all  three  kinetic  equations 
proposed  in  this  study,  Equations  3*12,  3*13  and  3*1*0  are 
shown  below. 

Kinetic  equation  AlL  =  H0(<rm-  a)(0  +  <r)r  (3.1*0 

d  t 


°  ( L, t) 


exp  td 


'm 


exp  td 


+ 


["  1  4-  Q'm 

L  4  J 


exp  ld  -  1 


r(L ,  7) 


_ exp  td 

exp  +  exp  I ^  —  1 


where  td  =  Horj  (  <rm  +  p)T 

L  =  H0LA<5  <rm0 

q 


(3-34) 


(3-35) 


Kinetic  equation 


(3.12) 
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*(L ,r) 


exp  td  -  1 _ 

exp  tjj  +  exp  l(-j  -  1 


(3-36) 


r(L,r) 


_ exp  td _ 

exp  +  exp  1^  -  1 


(3-37) 


where  td  =  Hor;T 

I.  =  Hp  A  0  q-mL 

°  q 

Kinetic  equation 


(3-13) 


a  (L,r) 


r(L,r) 


exp  td  -  1 


exp  tj  +  r 

1  -  tfmlexp  L  -  1 

L 

0  J 

_ exp  td 

exp  td  +  exp  ^  -  1 


(3-38) 

(3-39) 


where  td  =  H0  r j  (  0  -  ojr 

'd  =  HoLA0/3<rm  and  /3>(Tmalways 

q 

For  Kinetic  Equation  3*12,  only  simplification  1  is  required, 
while  for  Kinetic  Equation  3*13  "both  simplifications  1  and  2 
are  used.  Equations  of  this  type  have  been  derived  for  the 
flow  of  suspensions  through  porous  media  ( 15 ) » (44) , (45 ) , (46 ) . 

The  models  stated  above  are  derived  for  an  emulsion 
consisting  of  one  droplet  size,  flowing  in  a  linear, 
homogeneous  core  previously  saturated  with  water,  such  that 
all  the  droplets  in  the  emulsion  have  the  same  probability 
of  capture.  The  assumption  of  a  linear,  homogeneous  core  was 
required  so  that  <rm  could  be  assumed  constant  along  the 
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core,  in  order  to  carry  out  the  integration. 


3.4.5  N-Fractional  Size  Emulsion 


An  emulsion  is  generally  characterised  by  droplets  of 
various  sizes,  the  size  distribution  usually  being  described 
by  some  type  of  probability  distribution,  rather  than  of  one 
narrow  size  range.  Thus  the  probability  of  capture  for  each 
droplet  size  is  different  because  there  are  a  greater  number 
of  constriction  sites  available  within  the  core  for  the 
larger,  rather  than  smaller  droplets. 

The  above  models  for  the  single  sized  emulsion  can 
thus  be  extended  for  a  n-fractional  sized  emulsion  by 
introducing  the  variables  of  the  system  for  each  size  range. 
The  interaction  between  the  droplet  sizes  is  accounted  for, 
as  for  the  permanent  semi-permanent  capture,  by  modelling 
the  effective  capture  rate  for  each  droplet  size  range. 

Extension  of  the  model  for  Kinetic  Equation  3  •  14 
produces 


°  (L ,t) 


r(L,r) 


m ; 


exptdj  - 


1 


exp  tdj  + 


i  +  zmi 

.  (Jj  J 


exp  ld.  -  1 


exp  tdj 

exp  td  +  exp  ld.  -  1 


(3-40) 


(3.41) 


where  td  =  Hoj  rij  ( <rm-  +  0 ■ )  t- 

Id:  =  Hoj  A$0j  o-j  L 
J  P 
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and  subscript  j  refers  to  each  fractional  size  in  the 
emulsion . 

3.5  THE  EMULSION  FLOW  MODEL 

The  equations  reported  in  section  3*4.4  represent  the 
emulsion  flow  models  for  the  three  kinetic  equations  3*12, 
3.13  and  3*14  using  the  boundary  and  initial  conditions 
3.16,  3*17  and  3*27  and  simplifications  1  and  2. 

These  models  assume  that  the  rate  of  plugging  can  be 
modelled  by  a  rate  equation  and  the  rate  at  which  the 
retained  oil  saturation  increases  in  the  system  is  dependent 
on  the  saturation  within  the  system  and  the  concentration  of 
the  flowing  emulsion.  When  the  coefficients  of  the  model 
H  ,  (3  and  e  ,  have  been  evaluated  from  the  experimental 
data,  the  model  can  be  used  to  predict  the  flowing  emulsion 
concentration,  T  and  the  retained  oil  concentration,  o 
profiles  within  the  system  at  any  position  and  time. 

Some  of  the  properties  of  these  models  are: 

1)  the  models  allow  the  retained  oil  concentration  at 
the  inlet  end  of  the  system  to  increase  slowly 
with  time  to  its  maximum  concentration.  This 
boundary  condition  as  often  proven  difficult  to 
define  in  waterflooding  and  it  is  often  assumed  to 
jump  from  zero  to  the  maximum  value  at  the  instant 
flow  starts . 

2)  When  the  plugging  rate  constant  Hq  is  large,  the 
model  is  equivalent  to  the  moving  boundary  model. 
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3)  At  early  times,  the  model  allows  the  oil  front  to 
move  through  the  porous  medium  at  a  velocity 
greater  than  the  superficial  velocity,  while  at 
later  times,  the  front  can  fall  behind  the 
position  designated  by  the  superficial  velocity. 

4)  The  shape  of  the  front  is  controlled  by  the 
plugging  rate  constant  H  ,  while  the  position  of 
the  front  is  controlled  by  the  coefficients  6  and 

&  ■ 
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CHAPTER  4 

EXPERIMENTAL  EQUIPMENT  AND  PROCEDURE 

4.1  POROUS  MEDIA 

In  this  study  four  types  of  sandstone  porous 
media  were  used.  The  type  and  general  description  of  each  is 
given  below: 

1)  a  80  -  120  mesh  unconsolidated  Ottawa  sand  from 
Fisher  which  produced  cores  of  an  average  porosity  of  38%, 
with  a  permeability  range  of  14  -  18  darcies, 

2)  an  unconsolidated  Sundance  Kinsella  sand,  which 
produced  an  average  porosity  of  33% >  with  a  permeability 
range  of  155  -  320  millidarcies , 

3)  an  unconsolidated  Lloydminster  sand,  which 
depending  on  the  average  grain  size  produced  an  average 
porosity  of  37%  and  a  permeability  range  of  3*6  -  8.5 
darcies  and 

4)  a  consolidated  Berea  sandstone  with  a  porosity  of 
20%  and  a  permeability  of  39  millidarcies. 

Due  to  the  quantity  of  Lloydminster  sand  available, 
occasionally  the  sand  had  to  be  cleaned  and  reused.  The  sand 
was  initially  washed  with  toluene  to  remove  the  oil,  then 
oven  dried  at  120  °F  to  evaporate  the  toluene .  Finally  the 
sand  was  fired  at  500  °F  for  48  hours  to  remove  any 
organics  which  may  have  been  adsorbed  onto  the  sand. 
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4.2  EMULSION  AND  EMULSIFIER  SYSTEM 

The  emulsion  used  in  this  study  was  a  Triton  X-100 
surfactant  stabilised  emulsion  of  tapwater  and  Bon  Accord 
oil.  Triton  X-100  is  a  water  soluble  iso-octyl 
phenoxypolyethoxy  ethanol  linked  to  up  to  ten  molecular 
ethylene  oxide  units.  Its  molecular  weight  is  647  at  m  =  10. 
Triton  X-100  has  the  formula 


(  ch3  ) 


where  8  <  m  <  10 

and  was  supplied  by  Union  Carbide  and  manufactured  by  Rohm 
and  Haas . 

The  pertinent  properties  of  Bon  Accord  oil  are 
summarised  in  Table  4.1.  The  viscosities  were  measured  with 
a  reverse  flow  Cannon-Fenske  viscometer,  the  densities  with 
a  Haake  Density  meter  and  the  interfacial  tension  with  a 
spinning  drop  apparatus . 


TABLE  4.1 


PROPERTIES  OF  BON  ACCORD  OIL 


Temperature  (  °F  )  7 0  80  100  120  150 

Viscosity  (  mm^  s  ^  )  18.8  15*9  11.2  -  6.0 

Density  (  g/cc  )  0.881  0.878  0.870  O.863  0.852 


6.0 


Surfactant  Concentration  (  °/o  ) 


0  1 


Interfacial  Tension  (  dynes/  cm  )  9*0  0.36 


2t 

0.24 
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4.3  FLUID  INJECTION,  CORE  HOLDER  AND  PRESSURE  RECORDING 
SYSTEMS 

A  diagram  of  the  flow  system  is  shown  in  Figure  4.1. 
Central  to  the  system  was  a  Ruska  pump  capable  of  28  constant 
rates,  varying  from  2.5  -  55 0  cc/hr.  Two  steel  cylinders,  one 
two  litre  cylinder  containing  water  and  a  one  litre  cylinder 
containing  emulsion,  were  connected  to  the  pump.  The  pump 
injected  mercury  into  the  cylinder,  which  in  turn  displaced 
the  fluid  of  interest  into  the  core. 

The  core  holder  for  the  unconsolidated  sand  was  a  4.8 
cm.  I.D.  ,  30cm  long  PYC  tube  with  an  internal  length  of 
23.2  cm  between  the  end  block  flanges.  The  end  blocks  were 
made  of  PVC,  contained  400  mesh  screens  and  were  attached  to 
the  core  holder  by  bolts  and  sealed  with  'O'  rings.  The 
internal  volumes  of  the  end  blocks  were  known. 

The  size  of  the  unconsolidated  core  holder  was  chosen 
in  order  that  a  number  of  pore  volumes  of  emulsion  could  be 
injected  into  the  core  without  recharging  the  pump.  The  pore 
volume  of  the  core  holder  was  approximately  170  cc,  while 
the  maximum  throughout  of  the  pump  was  5°°cc  without 
changing  cylinders.  This  approach  was  carried  out  because  it 
was  felt,  and  subsequently  established,  that  stopping  and 
restarting  the  flow  would  have  a  disturbing  effect  on  the 
conditions  within  the  core. 

The  core  holder  for  the  unconsolidated  sand  was  a 
rubber  sleeve  enclosed  in  a  steel  vessel,  see  Figure  4.2. 

The  sleeve  was  squeezed  onto  the  core  by  pressuring  an  oil 


0-1000  psi  Heise  Gauge 
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FIGURE  4.1  SCHEMATIC  DIAGRAM  OF  THE  EXPERIMENTAL  EQUIPMENT 
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FIGURE  4.2  SCHEMATIC  DIAGRAM  OF  THE  CONSOLIDATED 

SAND  CORE  HOLDER 
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enclosed  within  the  vessel. 

The  pressure  drop  across  the  core  was  measured  just 
prior  to  the  inlet  to  the  core  holder.  Depending  on  the 
expected  pressure  drop,  either  a  1000  psi  Heise  gauge  or  a 
mercury  manometer  was  connected  to  the  flowline  to  provide 
a  direct  indication  of  the  pressure .  A  Validyne  differential 
pressure  transducer,  either  1,  5»  20  or  1000  psi  was  also 
connected  to  the  line  to  provide  a  continuous  recording  of 
the  pressure  history.  The  signals  from  the  transducer  were 
recorded  on  chart  paper. 

4.4  EMULSION  PREPARATION 

The  emulsion  was  prepared  in  400  cc  hatches.  Exactly 
10  cc  of  surfactant,  constituting  5%  of  the  emulsion 
prepared  in  one  hatch,  was  added  to  270  cc  of  tapwater  and 
120  cc  of  oil.  The  surfactant  concentration  used  in 
preparing  the  emulsion  was  well  in  excess  of  the  minimum 
amount  required  to  obtain  a  stable  emulsion.  The  whole 
mixture  was  subsequently  stirred  in  an  Osterizer  blender 
for  two  minutes.  The  produced  emulsion  was  then  diluted  down 
to  the  required  quality  by  adding  tapwater.  In  calculating 
the  emulsion  quality  (  r  ),  the  surfactant  Triton  X-100  was 
considered  to  be  part  of  the  aqueous  phase. 

4.5  PROPERTIES  OF  THE  EMULSION 
4.5*1  Density 


The  density  of  the  emulsion  was  measured  using  a 


>  ■ 
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Haake  Density  meter  and  was  0.981  g/cc  for  an  emulsion 
quality  of  0.15  and  0.963  g/cc  for  an  emulsion  quality  of 
0.30. 

4.5.2  Droplet  Size  Distribution 

The  droplet  size  distribution  of  the  emulsion  was 
measured  with  a  Coulter  Counter  Model  TAII,  using  a  130 
micron  aperature  tube .  As  the  measurement  of  oil-in-water 
macroemulsion  droplet  sizes  using  a  Coulter  Counter  is  now 
a  standard  procedure  and  found  in  a  number  of  papers  (12), 
(36)  and  (37) »  the  operation  of  this  device  will  not  be 
described.  It  will  however  be  noted  that  the  dilutent  used 
was  Isoton  and  the  procedure  used  was  that  described  in  the 
Coulter  Counter  Manual  (4?).  As  primary  interest  was 
placed  on  the  large  droplets  in  the  emulsion,  a  130  micron 
aperature  tube  was  chosen  so  that  the  maximum  droplet  size 
was  measureable  within  the  range  of  the  tube,  which  was  2.0 
-  40.3  microns . 

4.50  Emulsion  Stability 

Prior  to  the  emulsion  flow  study,  a  preliminary  study 
was  carried  out  to  find  a  stable  emulsion  by  changing  the 
parameters  of  oil  type,  surfactant  type,  concentration  and 
mixing  time .  During  this  experiment  it  was  found  that  the 
Zeta  meter,  a  device  for  establishing  emulsion  stability, 
gave  different  conclusions  to  those  from  the  degree  of 
separation  during  centrifution  of  the  emulsion.  It  was 
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therefore  assumed  that  for  some  reason,  the  Zeta  meter 
was  producing  erroneous  results.  Thus  no  quantitative 
estimate  of  emulsion  stability  was  obtained. 

For  the  emulsion  used  in  this  study,  no  oil  phase 
separation  occurred  over  a  period  of  weeks  when  stored  in 
an  air  tight  container.  However,  if  left  open  to  the 
atmosphere,  within  days,  oil  separated  out  from  the 
emulsion.  When  left  over  a  period  of  time  in  an  air  tight 
container,  the  emulsion  did  tend  to  cream,  but  it  was 
easily  redispersed  upon  shaking. 

4.5 .4  Viscosity 

The  viscosity  of  the  emulsion  was  measured  using  a 
reverse  flow  Cannon-Fenske  viscometer.  The  viscosity 
correlations  were  obtained  by  accurately  diluting  down  the 
original  emulsion  with  tapwater  to  produce  a  viscosity  oil 
concentration  profile,  producing  the  seven  data  points  for 
r  =  0.15,  0.12,  0.09,  0.075,  0.06,  0.03  and  0.0.  Using  the 
seven  data  points,  the  viscosity  correlation 

»e  =  "w  <  1  +  cir  +  c2r2  }  (2-1} 

was  fitted  to  the  data  on  the  computer,  by  evaluating  the 

best  values  for  the  constants  and  C^.  In  all  cases,  the 

error  in  the  correlation  was  smaller  than  the  degree  of 

scatter  within  the  experimental  evaluation  of  the  emulsion 

viscosity . 

The  viscosity  correlation  was  then  used  to  evaluate 
the  concentration  of  the  produced  emulsion  in  the 


■ 
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displacement  tests. 


4.6  TYPICAL  RUNS 

4.6.1  Core  Preparation 

Prior  to  each  run,  a  sieve  analysis  of  the  sand  was 
carried  out,  followed  by  recombining  the  sand.  The  cores 
were  packed  by  tamping  the  core  holders  with  a  hammer  until 
no  visible  settlement  occurred  while  obtaining  a  core  length 
of  23-4  cm.  The  mass  of  the  core  holder  was  determined 
before  and  after  it  was  packed. 

The  cores  were  then  evacuated,  saturated  with  tapwater 
and  reweighed.  From  the  series  of  weighings,  the  sand  weight, 
the  water  weight  and  thus  the  volume  of  water  within  the 
core  holder  were  determined.  From  the  volume  of  water  within 
the  core  holder,  the  pore  volume  and  porosity  of  the  core 
were  determined  by  subtracting  the  value  for  the  water 
enclosed  in  the  core  holder  end  blocks.  The  absolute 
permeability  was  determined  routinely  by  the  flow  of 
tapwater,  at  which  time  the  transducer  calibration  was 
checked . 

4.6.2  The  Emulsion  Flood 

Two  types  of  emulsion  flood  were  carried  out  in  this 
study.  Type  1  was  carried  out  to  investigate  objectives  1 
and  2  of  this  study  and  type  2  to  investigate  objectives  3 
and  4.  The  conditions  of  these  floods  are  described 
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separately  below: 

Type  1  : 

In  these  experiments,  the  emulsion  was  flooded 
at  a  constant  rate  through  the  core  until  steady  state 
conditions  were  achieved.  Steady  state  was  designated  by  no 
further  change  in  the  pressure  drop  across  the  core.  When 
steady  state  conditions  were  achieved,  the  pressure  drop, 
flow  rate  and  emulsion  viscosity  were  recorded,  from  which 
the  permeability  of  the  core  to  emulsion  was  determined. 

The  flow  rate  was  then  lowered  and  the  flooding 
continued  until  a  new  steady  state  condition  was  obtained. 
This  procedure  of  lowering  the  flow  rate  and  evaluating  the 
core  permeability  to  emulsion  at  steady  state  was  continued 
until  a  predetermined  lower  flow  rate  was  reached.  The 
procedure  was  then  reversed  and  the  permeability  was 
evaluated  through  a  series  of  increasing  flow  rates,  back 
to  the  original  flow  rate .  This  procedure  was  carried  out  to 
establish  that  the  core  properties,  such  as  wettability,  had 
not  changed  due  to  the  adsorption  of  organics  from  the  oil 
or  the  surfactant.  The  net  result  of  this  procedure  was  to 
determine  the  effect  of  emulsion  flow  rate  or  the  pressure 
gradient,  on  the  permeability  of  the  core  to  the  emulsion. 

The  lower  limit  of  the  flow  rate  was  determined  by  the 
accuracy  with  which  the  pressure  drop  could  be  evaluated 
using  the  transducers.  At  low  rates,  the  error  in  the 
measurement  was  large  in  comparison  to  the  magnitude  of  the 
pressure  drop  and  thus  could  give  a  misleading  estimate  of 


' 


59 


the  core  permeability  to  the  emulsion. 

Homogeneity  of  the  emulsion  throughout  the  run  was 
assured  by: 

1)  producing  all  the  emulsion  before  the  run  and 
storing  this  emulsion  in  one  storage  container  and 

2)  at  the  end  of  the  run  for  each  flow  rate,  the 
emulsion  cylinder  in  the  flow  system  was 
replenished  with  new  emulsion,  which  was  kept 
throughly  mixed  in  the  storage  container. 

This  technique,  rather  than  using  one  large  cylinder 
in  the  flow  system,  ensured  the  homogeneity  of  the  emulsion, 
because  in  a  large  cylinder  some  droplet  size  segregation 
does  occur  over  a  period  of  time.  Thus  the  droplets  would  be 
slightly  larger  than  average  at  the  beginning  of  the  run, 
becoming  progressively  smaller  with  time. 

In  these  experiments  the  produced  emulsion  was 
collected  in  one  large  container  until  steady  state  occurred, 
when  a  small  sample  was  collected  from  the  core  outlet  and 
used  to  evaluate  the  viscosity. 

Type  2  : 

In  this  experiment,  the  unsteady  state  time 
period  for  emulsion  flow  was  investigated.  After  the  absolute 
permeability  evaluation,  the  emulsion  was  injected  through 
the  auxiliary  lines  and  transducers  on  the  equipment  up  to 
the  valve  upstream  of  the  core  holder,  in  order  that  a 
material  balance  could  be  carried  out  on  the  core . 
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For  a  constant  injection  flow  rate,  the  produced  emulsion 
was  collected  in  5°  cc  tubes  which  were  placed  in  a 
fractional  collector.  The  viscosity  of  the  produced  emulsion 
was  then  evaluated  for  each  tube  and  in  one  experiment,  the 
droplet  size  distribution. 

From  the  emulsion  viscosity,  the  emulsion  concentration 
was  calculated,  from  which  a  material  balance  was  carried  out 
on  the  system  to  evaluate  the  retention  of  the  oil  within  the 
core .  The  results  were  then  converted  to  refer  to  the  ends 
of  the  porous  medium,  rather  than  the  inlet  and  outlet  of  the 
core  holder.  From  these  results,  the  outlet  emulsion 
concentration  history  and  the  oil  saturation  within  the  core 
were  evaluated. 

From  the  results  of: 

1)  maximum  oil  saturation  within  the  core, 

2)  produced  emulsion  concentration  history  and 

3)  flow  rate,  inlet  emulsion  concentration,  core 
porosity,  length  and  area 

the  mathematical  models  dsecribed  in  Chapter  3  were  fitted  to 
the  experimental  results  in  order  to  evaluate  the  constants 
H  ,  /3  and  e  .  A  computer  program  COMPLX  was  used  to  find  the 
best  values  for  these  constants.  The  computer  program 
minimises  the  error  in  the  fitting  of  the  model  to  the 
results  using  the  Simplex  method  of  Box  (48)  and  Nelder  and 
Mead  (49).  From  the  evaluated  constants  for  the  model,  the 
saturation  history  of  the  core  was  calculated  using  the 


model . 


CHAPTER  5 

RESULTS  AND  DISCUSSION 


5*1  SAND,  CORE  AND  EMULSION  PROPERTIES 

A  sieve  analysis  of  the  sands  used  in  this  study  is 
presented  in  Table  5*1*  The  core  properties  of  porosity  and 
permeability  for  each  run  are  presented  in  Table  5*2.  The 
emulsion  properties  of  viscosity  and  concentration  for  each 
run  are  also  shown  in  Table  5*2. 

Tables  5*3  and  5*4  show  the  particle  size  distribution 
data  as  obtained  from  the  Coulter  Counter  for  the  15  and  3° 
percent  concentration  emulsions,  respectively.  Examination 
of  these  results,  shows  that  channels  1  and  2  were  cut-off 
and  this  procedure  was  carried  out  to  overcome  background 
noise  which  effects  the  results  of  the  lower  channels. 
Therefore  channel  3  represents  the  droplet  size  range,  0  - 
2.52  microns  rather  than  the  range  2  -  2.52  microns.  Figures 
5.1  and  5*2,  represent  a  correlation  of  this  data  on 
semi-log  paper,  equivalent  to  the  paper  provided  by  Coulter 
Electronic  Inc.  (  form  4221030  ).  The  curves  have  the 
typical  form  of  the  bell  shaped  distribution  curve . 

Figure  5*3  shows  a  representation  of  the  cumulative 
frequency  curves  on  log-probability  paper.  A  straight  line 
on  this  type  of  paper  demonstrates  that  the  emulsion 
particle  diameter  is  log-normally  distributed.  The 
correlation,  Figure  5*3»  demonstrates  that  the  emulsions 
used  in  this  study  were  not  log-normally  distributed  but 
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RUN 

1 

100 

101 

102 

103 

104 

106 

109 

110 

111 

112 

113 

114 

115 

119 

120 

121 

123 

124 

125 

126 

12? 

128 

129 
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TABLE  5-1 

CORE  SAND  PROPERTIES 


WEIGHT  FRACTION  EACH  SIZE 
<45  >100  <100  >140  <140  >325  <325 


MESH 

MESH 

MESH 

MESH 

SAND 

SAND 

<354  >149  <149  >105  <105  >44 

<44 

WEIGHT 

TYPE 

MICRONS 

MICRONS 

MICRONS 

MICRONS 

GRAMS 

B 

SK 

- 

- 

- 

- 

716.7 

SS 

- 

- 

- 

- 

746.0 

SS 

- 

- 

- 

- 

746  .4 

SK 

- 

- 

- 

- 

742.4 

SS 

- 

- 

- 

- 

741 .4 

SK 

0.402 

0.220 

0.307 

0.071 

733.8 

LM 

0.?80 

0.120 

0.090 

0.010 

712.4 

LM 

- 

- 

- 

- 

691.8 

LM 

0  .780 

0.120 

0.090 

0.010 

711.2 

LM 

0.780 

0.120 

0.090 

0.010 

711.5 

LM 

0.780 

0.120 

0.090 

0.010 

711.2 

LM 

0.780 

0.120 

0.090 

0.010 

711.1 

LM 

0.780 

0.120 

0.090 

0.010 

710.4 

LM 

0.780 

0.120 

0.090 

0.010 

711.1 

LM 

0.780 

0.120 

0.090 

0.010 

711.3 

LM 

0.780 

0.120 

0.090 

0.010 

711.1 

LM 

0.780 

0.120 

0.090 

0.010 

711.3 

LM 

0.780 

0.120 

0.090 

0.010 

709.5 

LM 

O.562 

0.240 

0.180 

0.018 

711.1 

LM 

O.343 

O.36I 

0.28  5 

0.011 

710.5 

LM 

0.0 

0.651 

0.313 

0.038 

688  .8 

LM 

0.092 

0.408 

0.473 

0.027 

701.1 

LM 

0.497 

0.134 

0.348 

0.021 

717.8 
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TABLE  5-2 

CORE  AND  EMULSION  PROPERTIES 


Run 

Porosity 

Permeability 

darcies 

Emulsion 

Concentration 

Emulsion 

Viscosity 

mPa  .s 

Core 

Length 

cm 

Core 

Area 

2 

cm 

1 

0.205 

O.O39 

O.30 

2.15 

62.4 

20.2 

100 

0.359 

0.320 

0.30 

2.1 

23.4 

18.1 

101 

0.352 

18.9 

0.30 

2.3 

23.4 

18.1 

102 

0.346 

15.9 

0.30 

2.3 

23.4 

18.1 

103 

O.338 

0.155 

0.30 

2.1 

23.4 

18.1 

104 

O.343 

14.25 

0.30 

2.08 

23.^ 

18.1 

106 

0.333 

0.217 

0.30 

2.03 

23.4 

18.1 

109 

0.366 

6.45 

0.15 

1.46 

23.4 

18.1 

110 

0.368 

7.32 

0.15 

1.46 

23.4 

18.1 

111 

0.369 

7.64 

0.15 

1.50 

23.4 

18.1 

112 

0.376 

7.87 

0.15 

1.48 

23.4 

18.1 

113 

O.372 

8.11 

0.15 

1.49 

23.4 

18.1 

114 

0.375 

8.83 

0.15 

1.49 

23.4 

18.1 

115 

0.371 

8.06 

0.15 

1.50 

23.4 

18.1 

119 

0.377 

8.12 

0.15 

1.48 

23.4 

18.1 

120 

0.371 

9.50 

0.15 

1.50 

23.4 

18.1 

121 

0.352 

6  .50 

0.15 

1.51 

23.4 

18.1 

123 

0.379 

8.55 

0.15 

1.47 

23.4 

18.1 

124 

0.376 

8.22 

0.15 

1.46 

23.4 

18.1 

125 

0.375 

5.23 

0.15 

1.49 

23.4 

18.1 

126 

O.387 

5.46 

0.15 

1.49 

23.4 

18.1 

127 

0.399 

3.03 

0.15 

1.46 

23.4 

18.1 

128 

0.397 

3.64 

0.15 

1.51 

23.4 

18.1 

129 

0.379 

4.22 

0.15 

1.52 

23.^ 

18.1 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 
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TABLE  5-3 

DROPLET  SIZE  DISTRIBUTION 
15  PERCENT  EMULSION  CONCENTRATION 


MINIMUM  SIZE  DISTRIBUTION 

DROPLET  DIAMETER  PERCENT  GREATER  THAN 

Microns  DIFFERENTIAL  CUMULATIVE 


1.26 

- 

- 

1.59 

- 

- 

2.00 

18.8 

100.0 

2.52 

21 .4 

81.2 

3.17 

23.5 

59.8 

4.00 

19.6 

36.3 

5.04 

10.5 

16.7 

6.35 

2.6 

6.2 

0 

0 

0 

00 

0 . 6 

3.6 

10.08 

0.5 

3.0 

12.7 

0.5 

2.5 

16.0 

0.7 

2.0 

20.2 

0.7 

1.3 

25.4 

0.6 

0 . 6 

32.0 

0.0 

0.0 

40.3 

0.0 

0.0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 
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TABLE  5-4 

DROPLET  SIZE  DISTRIBUTION 
30  PERCENT  EMULSION  CONCENTRATION 


MINIMUM  SIZE  DISTRIBUTION 

DROPLET  DIAMETER  PERCENT  GREATER  THAN 

Microns  DIFFERENTIAL  CUMULATIVE 


1.26 

- 

- 

1.59 

- 

- 

2.00 

20.4 

100.0 

2.52 

22.5 

79.6 

3.17 

21.4 

57.1 

4.00 

14.4 

35.7 

5.04 

6  .6 

21.3 

6.35 

2.6 

14. 7 

8.00 

1.4 

12.1 

10.08 

3.0 

10.7 

12.7 

2.7 

7.7 

16 .0 

2.0 

5.0 

20.2 

2.0 

3.0 

25.4 

1.0 

1.0 

32.0 

0.0 

0.0 

40.3 

0.0 

0.0 

00 

90 

80 

70 

60 

50 

40 

30 

20 

10 

0 
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THE  15%  QUALITY  EMULSION 
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5.2  DROPLET  SIZE  DISTRIBUTION  FOR 
THE  30%  QUALITY  EMULSION 
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FIGURE  5.3  EMULSION  DROPLET  SIZE  DISTRIBUTION 
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had  a  higher  frequency  of  the  larger  droplets  than  that 
characterised  by  a  log-normal  distribution. 

From  Figure  5»3»  the  value  for  the  emulsion  median 
diameter  is  easily  identifiable .  The  median  diameter  is  the 
value  of  the  droplet  size  about  which  the  volume  of 
particles  in  the  distribution  is  divided  equally.  The  two 
emulsions  had  an  identical  median  diameter  of  3*^5  microns. 
Although  the  median  diameters  for  the  two  emulsions  were 
identical,  there  were  a  higher  proportion  of  larger  droplets 
in  the  higher  concentration  emulsion.  The  values  for  the  99 
percent  less  than  emulsion  droplet  size  were  estimated  from 
Figure  5*3  to  be  25 and  21  microns  for  the  3°  and  15 
percent  concentration  emulsions,  respectively.  For  both 
emulsions  there  were  no  droplets  larger  than  32  microns. 

5.2  STEADY  STATE  FLOW  CHARACTERISTICS  OF  EMULSIONS  IN  A 
POROUS  MEDIUM 

Twenty  runs  were  carried  out  to  investigate  the  steady 
state  flow  characteristics  of  emulsions  in  a  porous  medium. 
The  permeability  of  the  core  to  emulsion  is  shown  in  Table 
5.5  for  each  run.  The  fractional  permeability,  the 
permeability  of  the  core  to  emulsion  as  a  fraction  of  the 
absolute  permeability,  is  shown  in  Table  5-6.  A  number  of 
phenomena  relating  to  emulsion  flow  in  a  porous  medium  can 
be  observed  from  these  results. 


PERMEABILITY  OF  THE  CORE  TO  EMULSION 
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5.2.1  Reduction  of  the  Core  Permeability  by  Emulsion  Flow 


It  was  postulated  earlier  in  Section  2.2  of  this 
study,  that  the  flow  of  an  emulsion  through  a  porous  medium 
could  reduce  the  permeability  of  the  medium.  The  results 
reported  in  Tables  5*5  and  5-6,  show  that  this  was  true  for 
the  emulsions  used  in  this  study,  when  the  permeability  of 
the  core  was  less  than  9-5  darcies.  Two  parameters  which 
can  provide  an  estimate  of  whether  and  when  an  emulsion  can 
reduce  the  permeability  of  the  core,  are  discussed  later  in 
Section  5*3  of  this  study. 

5.2.2  Degree  of  Core  Permeability  Reduction  by  Emulsion  Flow 


McAuliffe  (12)  reported  that,  in  no  case  did  the 
emulsion  completely  plug  the  porous  materials  used  in  his 
study.  The  core  retained  one  to  ten  percent  of  the  initial 
permeability  and  the  emulsion  continued  to  flow  through  the 
core  . 

In  this  study  for  runs  1,  106,  100  and  103,  where  the 
core  permeabilities  ranged  from  39  millidarcies  to  320 
millidarcies ,  the  permeabilities  were  reduced  to  below  7 
millidarcies  in  all  cases.  These  values  represent  a 
fractional  permeability  of  less  than  5  percent  of  the 
initial  permeability.  For  all  these  runs,  steady  state  was 
never  attained  because  the  pressure  drop  across  the  cores 
approached  400  psi  for  runs  100,  103  and  106  and  1000  psi 


' 


for  run  1,  with  these  pressures  representing  the  limits  of 
the  respective  core  holders  used  in  the  runs.  The  number  of 
pore  volumes  flowed  through  the  cores  before  the  pressure 
limit  of  the  core  holders  were  reached  were  1.0,  1.5,  3*1 
and  1.5  for  runs  1,  100,  103  and  106,  respectively. 
Therefore,  because  the  steady  state  core  permeability  to  the 
emulsion  was  never  attained  and  the  core  permeabilities  were 
so  low  that  extremely  high  pressure  gradients,  hundreds  of 
psi/ft,  were  required  for  flow  at  reservoir  rates,  then  the 
condition  of  the  core  can  be  termed  as  completely  plugged. 

It  can  therefore  be  postulated  from  the  results  of 
this  study,  that  there  are  three  possible  degrees  of  core 
permeability  reduction  due  to  the  flow  of  an  emulsion 
through  a  porous  medium. 

Firstly,  and  characterised  in  runs  101,  103  and  104 
for  which  the  core  permeabilities  ranged  from  1 4.25  to  18. 9 
darcies,  the  emulsion  does  not  plug  the  porous  medium  and 
reduce  the  permeability.  In  these  runs,  the  permeability  of 
the  porous  medium  was  not  reduced  because  the  emulsion 
droplets  were  too  small  to  plug  the  pore  throats . 

Secondly,  and  characterised  in  runs  109  through  129 » 
for  which  the  core  permeabilities  ranged  from  3  to  9-5 
darcies,  the  emulsion  plugged  the  core,  partially  reducing 
the  permeability  of  the  core . 

Thirdly,  and  characterised  in  runs  1,  100,  103  and  106 
the  emulsion  completely  plugs  the  core,  reducing  the 
permeability  to  such  an  extent  that  there  was  essentially  no 
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flow  at  normal  reservoir  pressure  gradients. 

A  parameter  for  predicting  the  degree  of  permeability 
reduction,  whether  zero,  partial  or  complete,  is  discussed 
in  Section  5*3  of  this  study. 

Examination  of  the  results  shown  in  Table  5*5»  some  of 
which  are  shown  in  Figure  5*4,  demonstrates  that  in 
general,  for  an  unconsolidated  porous  medium,  a  core  with  an 
higher  initial  permeability  will  retain  a  higher 
permeability  after  emulsion  flow,  than  a  core  with  a  lower 
initial  permeability.  This  observation  is  not  true  in  every 
case  because  permeability  is  an  average  property  of  the 
porous  medium,  while  it  is  the  particular  combination  of  the 
pore  sizes  that  determines  the  permeability  of  the  core  to 
an  emulsion. 

Although  there  is  a  general  relationship  between  the 
permeability  of  the  core  to  emulsion  and  the  initial  core 
permeability,  the  results  which  are  shown  in  Figure  5*5 
demonstrate  that  there  is  no  simple  relationship  between  the 
degree  of  pore  plugging  or  permeability  reduction  and  the 
core  permeability.  A  comparison  of  the  data  for  runs  109, 

125  and  126  shown  in  Figure  5-5>  demonstrates  that  the 
fractional  permeability  of  the  core  to  emulsion  flow  was 
higher  for  the  lower  permeability  cores. 

The  results  of  runs  109,  125  and  126  demonstrate  that 
the  core  with  the  higher  initial  permeability  retained  a 
higher  permeability  after  emulsion  flow,  but  this 
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represented  a  lower  fractional  permeability  of  the  initial 
permeability.  Therefore  it  was  postulated  that  there  is 
no  simple  relationship  between  the  degree  of  permeability 
reduction  and  the  core  properties  of  porosity  and 
permeability.  There  is  no  simple  relationship  because 
porosity  and  permeability  are  volume  averaged  qualities  of 
the  pore  sizes,  while  it  is  the  particular  combination  of 
these  pore  sizes  that  determines  the  shape  of  the  fractional 
permeability  versus  pressure  drop  or  flow  rate  profile. 

The  smooth  nature  of  the  results  of  the  fractional 
and  actual  permeability  of  the  core  to  emulsion  versus 
flow  rate,  shown  in  Figures  5-^  and  5*5»  demonstrates  that 
no  changes  in  the  rock  properties  occurred  during  the  runs . 
These  results  imply  that  the  adsorption  of  the  surfactant 
or  oil,  if  they  occurred,  did  not  affect  the  plugging  of 
the  core . 

5.2.3  Non-Darcy  Flow  Nature  of  Emulsions  in  a  Porous  Medium 


In  section  2.2  of  this  study  it  was  postulated,  by 
examining  Equation  2.17,  that  emulsion  flow  is  a  non-Darcy 
flow  process  when  the  emulsion  can  plug  the  porous  medium. 
This  phenomena  was  postulated  because  emulsion  droplets  are 
deformable  and  for  a  higher  pressure  gradient  they  will 
deform  to  a  greater  extent  and  pass  through  a  smaller 
constriction.  Therefore  as  the  pressure  gradient  is  reduced, 
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the  pore  throat  which  can  be  plugged  will  increase  and 
the  permeability  of  the  core  to  emulsion  will  decrease. 

The  relationship  between  the  pressure  drop  and  the 
flow  rate  for  run  124,  which  is  a  representative  example  of 
a  core  which  was  plugged  by  emulsion  flow  in  this  study,  is 
shown  in  Figure  5*8.  The  circular  data  points  represent  the 
flow  of  water  through  the  porous  medium.  The  linear 
relationship  indicates  the  flow  of  water  obeyed  Darcy's  Law, 
with  the  permeability  of  the  core  being  8.22  darcies.  The 
dashed  line  on  Figure  5*8  represents  the  pressure  -  flow 
rate  relation  for  a  liquid  of  equal  viscosity  to  the 
emulsion  and  which  obeys  Darcy's  Law.  The  triangular  data 
points  represent  the  pressure  -  flow  rate  profile  for  the 
emulsion  used  in  the  run.  The  non-linearity  of  the  line 
through  the  triangular  data  points  indicates  that  emulsion 
flow  is  a  non-Darcy  flow  process  when  the  emulsion  plugs 
the  porous  medium.  The  results  for  this  run,  reported  in 
Table  5*8,  show  that  the  fractional  permeability  of  the 
core  changed  from  O.588  to  0.812  as  the  flow  rate  increased 
from  3°  "to  400  cc/hr. 

This  non-Darcy  flow  process  for  emulsions  was 
originally  demonstrated  by  McAuliffe  (12).  This  phenomenon 
was  not  observed  by  Alvarado  (34),  who  reported  that  the 
reduction  of  the  core  permeability  was  independent  of  the 
flow  rate  and  the  flow  of  an  emulsion  in  a  porous  medium 
was  classified  as  a  Darcy  flow  process. 

The  degree  of  this  non-Darcy  flow  process  for 
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FIGURE  5.6  EXAMINATION  OF  DARCY'S  LAW 
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emulsions  is  a  function  of  the  pore  size  distribution  of  the 
porous  medium  and  the  size  of  the  droplets  in  the  emulsion. 
Properties  conducive  to  observing  this  phenomenon  are : 

1)  a  sharply-distributed  pore  size  distribution  and 

2)  an  emulsion  containing  large  rather  than  small 
droplets  because  they  deform  to  a  greater  extent 
under  a  particular  pressure  gradient. 

The  results  of  runs  101,  102  and  104  demonstrate  that 
when  the  emulsion  does  not  plug  the  porous  medium,  then 
the  flow  of  an  emulsion  is  a  Darcy  flow  process.  The  non- 
Darcy  flow  charactersitic  of  emulsion  flow,  described  by 
McAuliffe  (12)  as  'pseudo  non-Newtonian',  is  therefore  only 
evident  when  the  emulsion  droplets  are  large  enough  to  plug 
the  porous  medium. 

5.3  PREDICTION  TECHNIQUES  FOR  THE  DEGREE  OF  CORE  PLUGGING 
d g/^e  Ratio 

Herzig  et  al  (15)  and  Sakthivadivel  (46)  worked 
on  the  process  of  deep  filtration,  in  which  suspensions  are 
mechanically  entrapped  within  a  porous  medium.  They  showed 
that  the  parameter  dg/ds»  th e  ratio  of  the  sand  grain 
diameter  to  the  maximum  particle  size  in  the  suspension, 
could  be  used  to  predict  whether  the  porous  medium  would 
become  plugged. 

Their  analysis  can  be  extended  to  predict  whether  an 
emulsion  can  plug  a  porous  medium  by  using  the  parameter 
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d  /d  ,  the  ratio  of  the  sand  grain  size  to  the  deformed  size 


of  the  largest  droplet  in  the  emulsion. 


the  ratio  of  the  sand  grain  size  to  the  maximum  oil  droplet 
size  in  the  emulsion  might  he  used,  however,  this  ratio 
would  not  he  truly  valid  as  emulsion  droplets  are 
deformable.  The  use  of  the  ratio  d ^/d^,  rather  than  d g/^e» 
would  provide  a  more  optimistic  estimate  of  whether  an 
emulsion  could  plug  the  porous  medium. 


The  ratio  d  /d  can  he  used  to  predict  whether  a  core 


will  he  plugged  hy  an  emulsion  hy  comparing  this  ratio  with 


the  ratio  of  the  sand  grain  size  in  the 


porous  medium  to  the  pore  throat  size. 

For  a  porous  medium  consisting  of  uniformly  sized 
spherical  particles,  the  ratio  of  the  grain  size  to  the 
minimum  pore  throat  and  maximum  pore  opening  in  the  porous 
medium,  for  the  various  packing  schemes  of  cubic,  hexagonal 
and  rhomhehedral  has  been  evaluated.  The  values  for  the 
various  packing  schemes  are  shown  in  Table  5*7* 


TABLE  5-7 


Minimum  Pore 
Throat  Diameter 


Maximum  Pore 
Opening  Diameter 


Packing 


Hexagonal 


Rhombic 


Cubic 


2.415 
6  .493 
6  .493 


1.366 
1 .894 
2.415 


. 
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Thus  for  spherical  grains  in  a  random  pack,  the 
minimum  pore  throat  size  is  given  by  0.154  (  1/6.493  )  of 
the  grain  size  and  this  represents  the  easiest  pore  which 
can  be  plugged.  Therefore  for  plugging  to  occur  in  a  random 
pack,  the  ratio  d g/de  has  to  be  smaller  than  6.493»  for  a 
single  droplet  to  plug  the  pore  constriction. 

In  a  random  pack,  the  maximum  pore  throat  size  is 
characterised  by  the  pore  throat  for  cubic  packing. 
Therefore,  for  complete  plugging,  the  ratio  d^/d^  has  to  be 
smaller  than  2.41 5*  Therefore  theoretically,  in  a  random 
core  pack  consisting  of  cubic,  hexagonal  and  rhombehedrally 
packed  spherical  grains,  then  for 

dg/dg  >  6.493  *  no  plugging 

2.415  <  dg/he  <  6.493  :  partial  plugging  (5*1) 

dg/dg  <  2.415  :  complete  plugging 

Herzig  et  al  (15)  stated  that  for  plugging  to  begin 

dg/de  <  6.493  :  plugging  with  one 

droplet  capture 

d/d  <  10.0  :  plugging  with  three 

S  e  (5.2) 

droplet  capture 

d ^/dg  <  12.2  :  plugging  with  four 

droplet  capture 

A  representation  of  their  theory  is  illustrated  in  Figure 

5-7. 


0.10 


0.154 


d  :  Diameter  of  Emulsion  Droplet 
dg  :  Diameter  of  Sand  Grain 


FIGURE  5.7  SIZE  OF  EMULSION  DROPLET 

RELATED  TO  TRIANGULAR 
CONSTRICTION  SITE  (15) 
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Sakthivadivel  (46)  in  his  work  stated  that  for 


d  /d 
g  e 

Vde 

5.0  < 

plugging  dependent 
on  the  porosity 


<  5*0  :  plugging  occurs 

>14.0  :  no  plugging 

<  14.0  :  plugging  or  no  (5-3) 


Calculation  of  d  : 

g 

The  value  for  d  ,  for  each  core,  was 

g 

evaluated  using  the  modified  Kozeny  equation 


180  (  1  -  </>  ) 


03 


K 


where  K  :  permeability  (  cm  ) 

d  :  grain  size  (  cm  ) 

g 

<t>  :  porosity 

which  is  valid  for  an  unconsolidated  porous  medium  made 
from  spherical  particles  of  one  distinct  size.  The  values 
of  d  evaluated  for  each  run  are  shown  in  Table  5*8. 


Calculation  of  d  : 

e 

In  suspension  flow,  the  value  for  the 
maximum  particle  size  is  constant,  however  because  emulsion 
droplets  are  deformable,  the  value  for  the  deformed  size  of 
the  largest  emulsion  droplet  must  be  evaluated  at  a  specific 
pressure  gradient.  In  this  study,  the  value  of  dg  and  the 
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TABLE  5.8 
PARAMETER  d  g/d0 


d^/dg  @  Plugging 

RUN  rQ  dg  ro/dg  0  15  Type 


microns 

microns 

psi/ ft 

psi/ft 

psi/ft 

101 

12.7 

180 

0.141 

0 

00 

7.20 

7.65 

NP 

102 

12.7 

171 

0.148 

6.73 

6 .84 

7.26 

NP 

104 

12.7 

164 

0.154 

6.48 

6  .58 

6.99 

NP 

120 

10.5 

118 

0.178 

5.62 

5.68 

5-92 

PP 

114 

10.5 

108 

0.194 

5.14 

5.20 

5.42 

PP 

123 

10.5 

104 

0.202 

4.95 

5.00 

5.22 

PP 

124 

10.5 

103 

0.204 

4.90 

4.96 

5.17 

PP 

119 

10.5 

106 

0.198 

5.05 

5.10 

5-32 

PP 

113 

10.5 

105 

0.200 

5.00 

5.05 

5.27 

PP 

112 

10.5 

101 

0.208 

4.81 

4.86 

5-07 

PP 

111 

10.5 

104 

0.202 

4.95 

5.00 

5.22 

PP 

110 

10.5 

102 

0.206 

4.86 

4.91 

5.12 

PP 

121 

10.5 

106 

0.198 

5.05 

5.10 

5.32 

PP 

109 

10.5 

97 

0.216 

4.62 

4.67 

4.86 

PP 

126 

10.5 

79 

0.266 

3.76 

3. 80 

3.96 

PP 

125 

10.5 

83 

0.253 

3.95 

3-99 

4.16 

PP 

129 

10.5 

65 

0.323 

3.10 

3.13 

3.26 

PP 

128 

10.5 

61 

0.344 

2.90 

2.94 

3.06 

PP 

127 

10.5 

59 

0.356 

2.81 

2.84 

2.96 

PP 

100 

12.7 

22 

1.154 

0.88 

O.90 

0.95 

PP 

106 

12.7 

22 

1.15^ 

0.848 

0.86 

0.92 

PP 

103 

12.7 

18 

1.411 

O.70 

0.71 

0.75 

PP 

86 


ratio  dg/dg  were  evaluated  for  the  pressure  gradients  of 

0  psi/ft,  1  psi/ft  and  5  psi/ft.  The  value  of  d  was 

0 

evaluated  using  the  equation 

dP  (  2 

—  L ,  =27  cos 0  <  — 

dL  d  (d 

e 

maximum  emulsion  droplet  radius 
wetting  angle,  assumed  to  be  90° 
interfacial  tension 

length  of  the  droplet  over  which  the  pressure 
gradient  applies 

pressure  gradient 

pressure  differential 

the  diameter  of  the  deformed  emulsion  droplet  of 
radius  r  ,  under  the  applied  pressure  gradient. 

In  order  to  evaluate  the  deformed  droplet  diameter,  dg 
the  length  over  which  the  pressure  gradient  applies  has  to  be 
be  estimated.  In  this  study,  the  pressure  gradient  was 
assumed  to  apply  over  the  total  length  of  the  deformed 
droplet.  The  deformed  droplet  length  was  evaluated  as  the 
addition  of  the  minimum  and  maximum  radii  of  the  droplet 
plus  the  length  of  a  conical  shape  between  the  end  radii. 
Figure  5*8  shows  a  schematic  diagram  of  the  deformed  emulsion 
emulsion  droplet.  The  calculation  methods  for  the  deformed 
droplet  length  and  diameter  and  the  value  of  d /d  are 


(5-5) 


nr 


where 


Ap  = 


m 


6 

7 

L 


dP 
■■  ■  • 

dL 

AP  : 
d_  : 
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=  2 


'm 


r2 

rm 


+  rp  +  rmrp 
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m 


+ 


m 


+  rr 


+  rmrp 


FIGURE  5.8  SCHEMATIC  DIAGRAM  OF  THE 

DEFORMED  EMULSION  DROPLET 
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reported  in  Appendix  A.  The  value  for  the  deformed  emulsion 
droplet  length  for  each  run  is  recorded  in  Table  5-9* 

The  value  of  the  maximum  emulsion  droplet  diameter 
used  in  the  calculation  of  the  parameter  d g/^e»  was  ‘the  99 
percent  less  than  size,  evaluated  in  Section  5*1*  This  value 
was  used  because  the  exact  maximum  droplet  size  could  not  be 
evaluated  from  the  Coulter  Counter  analysis. 


The  Parameter  d  /d  : 

g  e 


The  ratio  d^/d^  for  the  pressure 
gradients  0,  1  and  5  psi/ft  is  shown  in  Table  5*8.  From  the 
values  of  the  ratio  d^/d^  calculated,  the  type  of  plugging 
expected  was  determined  using  inequality  5* 1* 

The  values  of  the  ratio  d^/de  for  the  1  psi/ft 
pressure  gradient,  shows  that  for  runs  101,  102  and  104  the 
values  were  greater  than  6.493,  which  implies  the  porous 
medium  should  not  become  plugged.  For  runs  100,  103  and  106, 
the  ratio  d^/de  was  less  than  2.415,  which  implies  the 
porous  medium  should  be  completely  plugged.  While  for  runs 
109  -  129,  the  ratio  was  between  the  values  of  2.415  and 
6.493,  which  implies  the  porous  medium  should  be  partially 
plugged.  Therefore  the  degree  of  plugging  determined  from 
the  ratio  d^/dg  showed  a  perfect  correlation  with  the 
experimentally  observed  results. 

The  correlation  parameter  d^/dQ ,  therefore  yields  a 
simple  technique  for  estimating  whether  an  emulsion  will 
plug  an  unconsolidated  core  by  having  knowledge  of: 
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TABLE  5-9 

CRITICAL  CAPILLARY  NUMBER 


RUN 

PORE 

RADIUS 

rp 

DROPLET 

RADIUS 

r 

0 

DROPLET 

LENGTH 

L 

CRITICAL  CAPILLARY 

NUMBER 

101 

13.85 

12.7 

- 

102 

13.15 

12.7 

- 

104 

12.66 

12.7 

25.4 

8.02 

X 

-4 

10 

120 

9.08 

10.5 

22.4 

3.36 

X 

io~2 

114 

8.31 

10.5 

23.2 

5.03 

X 

10-2 

123 

8.01 

10.5 

23.5 

5.61 

X 

10-2 

124 

7.93 

10.5 

23.6 

5.65 

X 

10"2 

119 

8.16 

10.5 

23.3 

4.98 

X 

10~2 

113 

8.08 

10.5 

23.4 

5.24 

X 

10-2 

112 

7.77 

10.5 

23.7 

5.83 

X 

10-2 

111 

8.01 

10.5 

23-5 

5.15 

X 

10-2 

110 

7.85 

10.5 

23.6 

5.3^ 

X 

10"2 

121 

8.16 

10.5 

23.3 

4.27 

X 

10-2 

109 

7.47 

10.5 

24.0 

5.60 

X 

10-2 

126 

6.10 

10.5 

25.4 

7.53 

X 

10'2 

125 

6.38 

10.5 

25.1 

6.74 

X 

10"2 

129 

4.98 

10.5 

26.5 

8  .76 

X 

10-2 

128 

4.71 

10.5 

26.8 

7.91 

X 

10-2 

127 

4.54 

10.5 

27 .0 

6.95 

X 

io“2 

100 

1.73 

10.5 

36.4 

2.41 

X 

10"2 

106 

1.66 

10.5 

36.4 

1 .85 

X 

io“2 

103 

1.36 

10.5 

36.7 

1 .62 

X 

io-2 
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1)  the  core  permeability  and  porosity 

2)  maximum  oil  droplet  size  in  the  emulsion  and 

3)  the  interfacial  tension. 

Care  however  should  be  taken  around  the  inequality  values 
because  of  the  errors  and  assumptions  made  in  the  analysis: 

1)  the  estimation  of  the  maximum  droplet  size  in  the 
emulsion, 

2)  estimation  of  the  length  L^,  over  which  the 
pressure  gradient  applies  and 

3)  assuming  the  value  for  d  calculated  from  the 

§ 

permeability  and  porosity  using  Equation  5«^  is 
valid  for  the  core . 

Examination  of  the  values  for  the  ratio  dg/de  shown  in 
Table  5-8,  demonstrates  that  a  change  in  the  applied 
pressure  gradient  does  have  a  moderate  effect  on  the  ratio 
dg/dg.  This  change  in  the  deformed  emulsion  droplet  size  is 
the  basis  of  the  non-Darcy  flow  characteristic  of  emulsions 
described  earlier. 

Evaluation  of  the  parameter  d^/d^  i°r  the  work  of 
other  researchers,  shows  that  for  Uzoigwe's  work  (3 6)  on 
glass  bead  packs,  the  ratio  varied  from  35  to  126,  which 
implies  the  experimentally  observed  result  of  no  plugging. 
The  value  of  the  ratio  for  Cartmill's  work  (38)  is  2b  to  59  > 
which  predicts  no  plugging  should  occur.  However,  Cartmill 
observed  that  the  i  to  1-J  micron  sized  emulsion  droplets 
were  screened  out  by  the  37  to  88  micron  glass  beads.  The 
ratios  d  /de  and  d /d^  are,  however,  only  valid  for  plugging 
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due  to  capillary  forces,  while  Cartmill  suggested 
electrostatic  rather  than  capillary  forces  were  the  reason 
for  the  plugging. 

The  experimental  results  of  this  study  and  the 
parameter  d ^/dg ,  further  demonstrate  that  a  core  with  any 
permeability  can  be  plugged  as  long  as  the  emulsion  droplets 
are  large  enough.  McAuliffe  (12)  stated  that  emulsions 
could  effectively  reduce  the  permeability  of  the  core  when 
the  initial  permeability  was  two  darcies  or  less. 


Capillary  Number 


The  capillary  number,  N  is  most  often 

c 


defined  as 


v 


N 


7 


(5.6) 


which  on  combining  with  equations  2.1 7  and  2.18  can  be 
written  as 


q  fi 


N  = 

c  A<t>y 


K  K  dP 
r 

7  <j)  dL 


K  K  2  cos 6  (  1  1 

r 


0  L, 


(5-7) 


r  r 

p  O' 


The  values  of  the  capillary  number  at  which  plugging  of  the 
porous  medium  commences  can  be  calculated  using  the  equation 

1 


N  = 
c 


K  K  2  cos  0  (  1 
r 


0  L 


(5-8) 


r  r 
P  o' 


In  this  analysis,  the  value  of  the  length  over  which  the 
pressure  gradient  applies  was  calculated  as  before,  the 
relative  permeability  is  taken  to  be  one,  because  there  is 
no  plugging  and  the  contact  angle  is  assumed  equal  to  9° 
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degrees.  Equation  5*8  can  therefore  he  written 


2  K 


N 


</>  L 


1 

r 


1 

r 


(5-9) 


The  value  for  the  pore  throat  size  was  taken  as  0.077  d  and 

o 

the  value  for  the  maximum  droplet  size  was  again  taken  as 
the  99  percent  less  than  size.  The  values  olf  the  capillary 
number  calculated  using  this  equation  represent  the  value  at 
which  plugging  of  the  porous  medium  will  commence .  When 
the  capillary  number  exhibited  in  the  flow  system  is  less 
than  this  critical  capillary  number,  then  the  oil  droplets 
in  the  emulsion  will  plug  the  porous  medium. 

The  values  of  the  critical  capillary  number  calculated 
using  Equation  5*9  are  reported  in  Table  5*9*  A 
representation  of  the  critical  capillary  number  versus  the 
core  permeability,  with  a  cross  plot  of  the  core  porosity, 
calculated  using  the  analysis  of  Equation  5*9»  is  shown  in 
Figure  5«9*  The  results  calculated  from  this  simple  model 
study,  show  that  the  critical  capillary  number  necessary 
for  pore  plugging  is  lower  for  the  low  and  high  permeability 
cores  and  thus  pore  plugging  will  commence  at  a  lower  flow 
rate .  Although  the  critical  capillary  number  is  lower  for 
the  low  and  high  permeability  cores,  the  pressure  drop 
necessary  for  no  plugging,  increases  as  the  permeability 
decreases.  The  critical  pressure  drop  and  flow  rate 
necessary  for  no  plugging  to  occur,  are  reported  in  Table 
5.10.  The  capillary  number  concept  might  therefore  be  a 
useful  parameter  for  determining  whether  an  emulsion  will 
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FIGURE  5.9  CRITICAL  CAPILLARY  NUMBER  FOR  PERMEABILITY  REDUCTION 


RUN 

101 

102 

104 

120 

114 

123 

124 

119 

113 

112 

111 

110 

121 

109 

126 

125 

129 

128 

127 

100 

106 

103 
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TAELS  5-10 


CRITICAL  BLOCKAGE  PARAMETERS 


CRITICAL  CAPILLARY  CRITICAL  CRITICAL  PRESSURE 

NUMBER  FLOW  RATE  VELOCITY  GRAT'IEUT 

cc/hr  ft/day  psi/ft 


8.02 

X 

10' 

-4 

207 

26.2 

0 

.2 

3.56 

X 

10' 

-2 

1 

•  3 

X 

104 

1530 

14 

.1 

5.03 

X 

10' 

-2 

1 , 

.98 

X 

104 

2290 

22 

.9 

5.61 

X 

10' 

-2 

2, 

.26 

X 

4 

10 

2600 

26 

.7 

5.65 

X 

10' 

-2 

2 , 

.28 

X 

4 

10 

2640 

27 

.8 

4.98 

X 

10' 

-2 

1 , 

.99 

X 

4 

10* 

2290 

24 

.8 

5.24 

X 

10' 

-2 

2, 

■  05 

X 

4 

10 

2390 

25 

0 

.8 

5.83 

X 

10' 

-2 

2 , 

.32 

X 

4 

10* 

26SO 

29 

•  9 

5.15 

X 

10' 

-2 

1 

■  98 

X 

4 

10* 

2340 

26 

.7 

5-3^ 

X 

10' 

-2 

2 

.10 

X 

4 

10* 

2490 

28 

.8 

4.2? 

X 

10' 

-2 

1 , 

.55 

X 

4 

10 

1910 

24 

.8 

5.60 

X 

10' 

-2 

2, 

.19 

X 

4 

10* 

2610 

34 

.1 

7.53 

X 

10' 

-2 

3. 

■  05 

X 

4 

10* 

3430 

57 

.4 

6.74 

X 

10' 

-2 

2, 

.65 

X 

4 

10* 

30s  0 

51 

.9 

8.76 

X 

10' 

-2 

3, 

•  51 

X 

104 

4030 

84 

•  5 

7.91 

X 

10' 

-2 

3. 

.25 

X 

1c4 

3560 

92 

.7 

6 . 96 

X 

10' 

-2 

2, 

.97 

X 

4 

10 

3230 

98 

.4 

2.4l 

X 

10' 

-2 

6, 

.46 

X 

4 

10* 

7830 

291 

1.85 

X 

10' 

-2 

4, 

.74 

X 

4 

10* 

6190 

305 

1 .62 

X 

10' 

-2 

4, 

.07 

X 

4 

10* 

5240 

379 
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plug  a  porous  medium  at  a  particular  flow  rate . 

The  values  for  the  critical  capillary  number 

_2 

calculated  in  this  study  varied  from  1.81  x  10  to  8.24  x 
_2 

10  for  the  cores  which  were  plugged.  Morrow  (32)  in  his 
study  on  the  mobilisation  of  residual  oil  reported 
capillary  numbers  of  1.43  x  10-2  to  1.68  x  10-2  for  similar 
oil  droplet  to  grain  size  ratios. 

For  runs  101  and  102,  the  capillary  number  was  not 
calculated  because  the  emulsion  droplet  was  smaller  than 
the  pore  throat  size  and  therefore  the  core  could  not  be 

plugged.  For  run  104,  the  capillary  number  was  evaluated  to 

-4 

be  7*02  x  10  .  This  capillary  number  predicts  that  the  core 

would  be  plugged  for  flow  rates  less  than  180  cc/hr.  Within 
the  experimental  accuracy  of  measuring  the  pressure  drop 
across  the  core,  there  was  no  plugging  of  the  core  for  the 
flow  rates  of  560,  320  and  100  cc/hr.  However,  when  the  pore 
throat  radius  and  the  emulsion  droplet  radius  are  similar, 
the  evaluation  of  the  critical  capillary  number  and  thus  the 
critical  flow  rate,  are  very  sensitive  to  the  input  data.  A 
sensitivity  analysis  on  the  critical  capillary  number  is 
reported  in  Appendix  B,  along  with  the  conditions  when  the 
capillary  number  concept  is  reliable  for  predicting  the 
plugging  of  a  core . 

For  runs  109  to  129,  the  capillary  number  during  the 

-  S  -3 

displacement  tests  varied  from  55  x  10  J  to  1.4  x  10  , 

while  the  capillary  numbers  varied  from  1.81  x  10  to  8.24 

_  p 

x  10  .  Therefore  for  the  runs  109  to  129,  the  capillary 


. 
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number  was  lower  than  the  critical,  verifying  that  plugging 
should  occur. 

Thus  although  there  are  errors  and  assumptions  made 
in  evaluating  the  capillary  number  and  measuring  the 
permeability  of  the  core  to  emulsion,  the  concept  of  using 
the  critical  capillary  number  for  predicting  whether  a  core 
will  be  plugged  by  an  emulsion  would  appear  to  be  a  useful 
concept . 

5.4  EMULSION  VISCOSITY  CORRELATION 

In  this  study,  the  emulsion  viscosity  correlation 

2 

<4  =  kw  (  i  +  c1r  +  c2r  )  (2.i) 

was  used  to  evaluate  the  produced  emulsion  concentration. 
From  the  produced  emulsion  concentration,  a  material  balance 
was  carried  out,  from  which  the  oil  retention  in  the  core 
was  established.  The  use  of  a  viscosity  correlation  for 
establishing  a  material  balance  was  used  by  Holm  (5°),  in 
the  investigation  of  the  retention  of  polmers  during  a 
flood . 

The  values  of  the  constants  C^  and  C^  for  the  runs 
119>  120,  128  and  129»  for  which  the  correlation  was  used, 
are  reported  in  Table  5 .11.  A  representative  figure  of  the 
correlation  versus  the  experimental  data  is  shown  in  Figure 
5*10  for  run  119-  In  all  cases,  the  error  in  the  correlation 
was  less  than  the  error  associated  with  measuring  the 
emulsion  viscosity. 
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FIGURE  5.10  EMULSION  VISCOSITY  CORRELATION 


RUN  119 
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TABLE  5-11 

CONSTANTS  FOR  THE  VISCOSITY  CORRELATION 
Viscosity  Constants 


Run 

mPa .  s 

2  -1 

mm  s 

A 

C2 

119 

1 .48 

1 .50 

2.738 

4.444 

120 

1  .50 

1.53 

2.967 

3-556 

128 

1  .51 

1 .54 

2.920 

4.622 

129 

1  .52 

1.55 

2.843 

5 .649 

It  was  therefore  believed  that  the  use  of  the 
viscosity  correlation  to  calculate  the  emulsion 
concentration  was  more  accurate  than  the  other  methods 
which  might  have  been  used: 

1)  distillation  of  the  produced  emulsion, 

2)  use  of  the  Coulter  Counter  and 

3)  the  electrical  resistivity  of  the  solution. 


5-5  UNSTEADY  STATE  FLOW  CHARACTERISTICS  OF  EMULSIONS  IN  A 
POROUS  MEDIUM 

Four  runs  119,  120,  128  and  129  were  carried  out  to 
investigate  the  unsteady  state  flow  characteristics  of 
emulsion  flow  through  the  porous  medium.  The  raw  and 
processed  data  used  to  calculate  the  produced  emulsion 
concentration  profiles  for  the  four  runs,  is  reported  in 
Appendix  C .  A  summary  of  some  of  the  pertinent  data 
collected  in  these  runs  is  shown  in  Tables  5*12  and  5* 13* 
The  produced  emulsion  concentration  profiles  for  the  runs 
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are  shown  in  Figures  5*11  to  5«1^« 

Explanations  based  upon  the  current  emulsion  and 
immiscible  flow  theories  are  proposed  to  describe  the 
findings  of  this  study.  The  interpretations  are  however 
limited  to  general  concepts  due  to  the  wide  differences  in 
the  properties  of  the  porous  media  and  the  flow  rates  used 
in  the  runs. 

5*5 -1  Produced  Emulsion  Concentration  Profile 

Although  it  was  not  possible  to  ascertain  the  exact 
breakthrough  time  for  each  run  in  this  study,  it  was 
observed  that  breakthrough  occurred  earlier  in  those 
experiments  conducted  at  lower  flow  rates.  An  exact 
breakthrough  time  could  not  be  evaluated  because  the 
produced  fluids  passed  from  clear,  to  cloudy,  to  a  distinct 
colouration  over  a  period  of  a  few  cetilitres.  For  all  the 
runs,  the  oil  phase  breakthrough  occurred  before  one  pore 
volume  of  emulsion  was  injected.  It  is  believed  that  one 
criterion  for  the  oil  phase  breakthrough  time,  is  the 
relationship  between  the  maximum  oil  phase  saturation  in 
the  core  and  the  injected  emulsion  concentration. 

In  the  low  concentration  emulsion  floods  of  Devereux 
(39)  the  maximum  oil  saturation  in  the  porous  medium  was 
many  times  the  inlet  concentration.  In  these  runs,  initially 
only  the  original  water  phase  in  the  core  and  the  external 
water  phase  of  the  emulsion  were  produced.  The  oil  phase  of 
the  emulsion  did  not  breakthrough  until  a  number  of  pore 
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FIGURE  5.11  PRODUCED  EMULSION  CONCENTRATION  :  RUN  119 
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FIGURE  5.12  PRODUCED  EMULSION  CONCENTRATION  :  RUN  120 
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FIGURE  5.13  PRODUCED  EMULSION  CONCENTRATION  :  RUN  128 
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FIGURE  5.14  PRODUCED  EMULSION  CONCENTRATION  :  RUN129 
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volumes  of  emulsion  were  injected.  A  possible  explanation 
for  the  delayed  oil  phase  breakthrough  is  the  high 
efficiency  of  the  porous  medium  at  screening  out  the 
emulsion  droplets.  In  this  study,  see  Table  5* 12,  the  oil 
saturation  SQ  was  only  marginally  greater  than  the  inlet 
concentration  of  the  emulsion  and  the  lower  screening 
efficiency  of  the  porous  medium  led  to  an  earlier  oil 
phase  breakthrough. 

It  is  believed  that  another  criterion  for  the  number 
of  pore  volumes  injected  before  oil  phase  breakthrough 
occurs  in  the  laboratory,  is  the  length  of  the  core.  This 
criterion  will  be  discussed  further  in  Section  5*6.3  of 
this  study. 

The  results  of  this  study  also  showed  that  a  steeper 
front  occurred  in  the  higher  flow  rate  runs.  This  can  be 
observed  by  comparing  the  produced  emulsion  profiles 
shown  in  Figures  5*11  to  5*1^* 

5.5.2  Oil  Retention  in  the  Porous  Medium 

A  comparison  of  the  magnitudes  of  the  maximum  oil 
saturation  within  the  cores  for  each  run  is  not  justifiable, 
as  the  properties  of  the  porous  medium  and  the  flow  rates 
varied  for  each  run.  However,  the  results  do  show  that  the 
oil  accumulation  in  the  core  for  all  the  runs  was  slightly 
greater  than  the  inlet  concentration,  see  Table  5*12.  Thus 
the  blockage  coefficient  (b)  for  each  run  was  always  less 
than  i/r,  which  for  these  runs  was  6.667. 
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The  results  of  this  study  also  confirm  the  observation 
of  McAuliffe  (12)  that  emulsified  oil  is  more  efficient  at 
reducing  the  permeability  than  oil  left  after  waterflooding. 
The  results  for  run  129*  reported  in  Appendix  C,  show  that  a 
plugged  oil  saturation  of  12  percent  reduced  the  core 
permeability  by  45.6  percent,  while  for  residual  oil  it 
would  be  approximately  25  percent  using  McAuliffe ' s  analysis 
(12)  . 

In  order  to  investigate  the  plugging  capability  of 
the  emulsion,  a  droplet  size  distribution  was  carried  out  on 
the  produced  emulsion  in  run  129*  The  results  of  this 
analysis  are  reported  in  Appendix  D,  and  a  smoothed  graph  of 
the  droplet  size  versus  pore  volumes  injected  is  shown  in 
Figure  5- 15-  The  size  distribution  of  the  emulsion  was 
divided  up  into  four  sizes  in  order  to  simplify  the 
manipulation  and  presentation  of  the  data. 

Comparison  of  the  produced  emulsion  droplet  size 
distribution  near  the  end  of  the  run  and  the  inlet  size 
distribution,  shows  that  the  characteristics  of  the  produced 
emulsion  were  approaching  the  inlet  conditions.  This  trend 
of  the  produced  emulsion  droplet  size  towards  the  inlet 
size  distribution,  implies  that  for  steady  state  conditions, 
the  conservation  of  droplet  size  within  the  emulsion  would 
occur.  This  phenomenon  of  no  aggregation  or  breaking  up  of 
the  droplets  at  steady  state  flow  conditions  in  a  porous 
medium,  for  low  concentration  emulsions,  was  also  observed 
by  Uzoigwe  (36 ) . 
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FIGURE  5.15  DROPLET  SIZE  DISTRIBUTION  OF  PRODUCED  EMULSION 


The  results  reported  in  Table  5-13  show  that  there 
were  different  degrees  of  retention  in  the  core  for  the 
different  droplet  size  ranges.  Within  the  limits  of  the 
evaluation,  the  oil  droplets  represented  in  channels  6  and 
7  +  of  the  Coulter  Counter  analysis  had  a  higher 
concentration  in  the  core  than  their  inlet  concentration. 
Channels  3+4  and  5  were  both  retained  with  a  concentration 
close  to  their  inlet  concentration. 

One  explanation  for  the  different  retentions,  for 
each  droplet  size  range,  is  that  the  larger  droplets  have 
more  plugging  sites  to  occupy  than  the  smaller  droplets. 


described  earlier  in 


section  5 • 3 >  then  the  droplets  of  channels  3+4,  5  and  6 
would  not  normally  be  retained  in  the  core  if  they  were 
the  only  droplets  flowing  in  the  core .  The  value  of  the 
ratio  dg/d^  for  each  channel  is  shown  in  Table  D4 .  Thus  it 
is  postulated  that  the  larger  droplets  plug  off  the 
constrictions  and  the  smaller  droplets  which  would  normally 
not  be  retained,  are  retained  in  cavern  sites  behind  the 
larger  droplets  or  in  combination  with  the  larger  droplets 
to  plug  off  a  constriction. 

Although  the  smaller  droplets  reached  their  inlet 
concentration  rapidly,  it  is  believed  that  they  were  still 
being  trapped  in  constrictions  along  with  larger  droplets 
throughout  the  run.  However,  the  droplets  were  breaking 
away  from  other  sites  and  because  the  droplets  were 
retained  at  a  concentration  similar  to  their  inlet 
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concentration,  their  concentration  in  the  core  remained 
approximately  constant. 

Examination  of  the  results  for  channels  7  and  8  near 
the  end  of  the  run,  see  Table  D2,  shows  that  these  channels 
had  a  larger  concentration  of  droplets  leaving  the  core 
than  were  being  injected.  A  possible  explanation  for  why 
these  droplets  were  leaving  the  core  at  a  higher  rate  than 
their  injection  rate,  is  that  the  pressure  drop  had  increased 
sufficiently  to  squeeze  the  droplets  through  the 
constrictions  and  because  the  droplets  were  retained  within 
the  core  at  a  concentration  higher  than  their  inlet 
concentration,  they  were  produced  at  a  higher  concentration. 

Thus  for  droplets  within  a  certain  size  range,  such  as 
was  exhibited  by  channels  7  and  8,  their  concentration  at 
steady  state  conditions  may  not  represent  their  highest 
concentration  recorded  throughout  the  run.  Their  maximum 
concentration  is  reached  prior  to  steady  state  conditions 
and  the  concentration  then  becomes  smaller  as  the  droplets 
leave  their  constriction  sites.  Thus  for  each  droplet  size 
there  is  an  equilibrium  concentration  at  steady  state 
conditions  and  a  maximum  concentration  and  these  conditions 
do  not  have  to  be  equal.  Thus  superimposed  on  the  rate  at 
which  droplets  are  plugged,  is  the  rate  at  which  the 
droplets  are  released  from  their  capture  sites. 

Thus  the  results  of  run  129,  shown  in  Tables  D2  and 
5.13  and  in  Figure  5-1 5,  show  that: 

l)  each  droplet  size  takes  a  different  time  for  the 
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produced  concentration  to  reach  the  inlet 
concentration, 

2)  there  are  different  degrees  of  retention  for  each 
droplet  size  range, 

3)  the  maximum  concentration  for  a  particular  droplet 
size  does  not  have  to  be  equal  to  the  concentration 
at  steady  state  conditions  and 

4)  the  produced  concentration  for  a  particular  droplet 
size  can  be  greater  than  the  inlet  concentration. 

5*5*3  Pressure  Profiles 

The  pressure  profiles  for  the  runs  119»  120,  128  and 
129  are  shown  in  Figure  5*16.  These  profiles  show  two 
characteristics  of  liquid-liquid  displacements  in  a  porous 
medium : 

1)  the  curved  pressure  profile  after  breakthrough  and 
before  steady  state  conditions  are  reached  and 

2)  the  curved  nature  of  the  pressure  drop  profile 
versus  time,  prior  to  breakthrough,  which 
designates  that  a  stabilised  front  had  not  formed. 

5*5*4  Effect  of  Interrupting  the  Flow  on  the  Equilibrium 

of  the  Flow  System 

Three  parameters,  the  pressure  drop,  the  produced 
emulsion  concentration  and  the  droplet  size  distribution 
were  examined  to  establish  the  effect  of  interrupting  the 
emulsion  flow  on  the  system  equilibrium. 
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The  pressure  profiles  shown  in  Figure  5*16, 
demonstrate  that  after  interrupting  the  flow,  it  took 
approximately  one  third  of  a  pore  volume  (  50  cc .  )  in  the 
small  core  holders,  to  reattain  the  pressure  drop  present 
before  interrupting  the  flow.  This  phenomenon  could  be  due 
to  the  emulsion  droplets  breaking  away  from  their  plugging 
sites  when  the  flow  was  stopped  and  on  restarting  the  flow, 
some  of  the  droplets  flowed  out  of  the  area  by  a  route 
other  than  their  plugging  site.  It  then  takes  some  time 
for  the  droplet  to  find  another  plugging  site. 

A  similar  disturbance  was  not  observed  in  the  emulsion 
viscosity  or  droplet  size  distribution  data,  although  this 
could  have  been  masked  due  to  the  influence  of  the  end 
blocks  and  the  size  of  the  samples  taken. 

5.6  EMULSION  FLOW  MODEL 

5*6.1  Emulsion  Flow  Model 

The  curves  on  Figure  5*17  to  5*20  represent  the  best 
fit  of  emulsion  flow  model  3,  Equation  5*1^»  to  the 
experimental  data  for  runs  119,  120,  128  and  129, 
respectively.  The  best  fit  values  for  the  coefficients  H  , 

(3  and  e  for  the  model  are  shown  in  Table  5«1^»  The  results 
demonstrated  that  the  model  can  effectively  describe  the 
early  time  periods  of  the  experimental  curves.  However,  the 
model  did  predict  an  earlier  time  to  reach  steady  state 


conditions . 
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FIGURE  5.17  EXAMINATION  OF  MODEL  3  :  RUN  119 
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FIGURE  5.18  EXAMINATION  OF  MODEL  3  :  RUN  120 
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FIGURE  5.19  EXAMINATION  OF  MODEL  3  :  RUN  128 
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FIGURE  5.20  EXAMINATION  OF  MODEL  3  :  RUN  129 
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TABLE  5.14 


BEST  FIT  COEFFICIENTS  FOR  MODEL  3 


RUN 

H 

0 

-1 

sec 

0 

R>t 

sec 

e 

a 

cm~l 

119 

0.812 

0.0764 

746 

0.727 

35.6 

120 

0.422 

0.0979 

753 

0.745 

18 .2 

128 

0.033 

0.136 

3859 

0.638 

8.5 

129 

0.029 

0.116 

3807 

0.659 

7-3 

It  was  found  in  this  study  that  model  3  fitted  the 
experimental  data  the  Lest,  followed  hy  model  1,  then  2  and 
because  of  this,  the  fits  of  models  1  and  2  are  not  shown. 
An  insight  into  why  model  3  fitted  the  data  the  best  can  be 
obtained  by  examining  the  coefficient  (3 . 

When  comparing  the  fits  of  models  1  and  2  together, 
the  value  of  (3  in  model  2  always  increased  to  the  upper 
limiting  value  set  in  the  computer  program.  In  model  2,  as 
(3  becomes  larger,  the  form  of  model  2  approaches  model  1, 
because  the  variable 

(  (3  -  o  )  approaches  (3 
and  the  function 

const.  (  am  -  o  )  (  (3  -  o  )  approaches  const.  (  crm  -  o  ) 
which  is  model  1. 

Comparing  models  3  and  1,  the  values  of  (3  in  model  3 
were  always  similar  to  a  .  If  model  1  was  the  better  model, 
then  (3  in  model  3  would  try  to  be  as  large  as  possible  in 
order  that  the  form  of  model  3  would  approach  model  1. 
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However,  (3  was  similar  to  <7m  ,  and  as  such  the  form  of  the 
rate  equation  is  convex  in  shape,  rather  than  linear.  Figure 
3-2  shows  the  comparative  shapes  of  the  different  kinetic 
rate  equations . 

The  method  by  which  the  model  is  fitted  to  the 
experimental  data  is  described  in  Appendix  E.  In  evaluating 
the  best  fit  coefficients  of  the  model,  it  was  found  that 
there  was  not  one  unique  set  of  constants  which  produced 
the  best  fit  to  the  data.  Due  to  the  interelation  between 
H  and  (3  in  the  model,  there  were  a  series  of  similar  values 
of  Hq  and  (3  which  provided  an  identical  fit.  Although  these 
coefficients  gave  the  same  fit  to  the  experimental  data,  it 
was  found  that  for  later  times,  some  of  these  sets  of 
coefficients  seriously  violated  the  material  balance  for 
the  system.  Thus  for  further  work  in  this  area,  longer  cores 
should  be  used  because  the  evaluation  of  the  coefficients 
of  the  model  is  more  sensitive  on  a  fully  developed  front. 

The  rate  coefficient  a  is  defined  as  HQA</>/q  and 
Trzaska  (45),  in  his  work  on  filtration,  stated  that  an 
increase  in  the  flow  rate  resulted  in  an  increase  in  the 
value  of  the  rate  coefficient  a  .  This  observation  states 
that  for  a  higher  flow  rate,  the  rate  of  plugging  increases 
and  a  steeper  front  exists  within  the  core.  Although  the 
core  properties  of  each  run  were  different  in  this  study, 
it  can  be  seen  in  Table  5*14  that  the  value  of  the  rate 
coefficients  a  were  higher  for  the  higher  flow  rate  runs. 

For  run  128,  the  model  did  not  provide  a  good  fit  to 
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the  experimental  data.  One  possible  reason  for  the  poor  fit, 
is  that  emulsions  are  made  up  of  droplets  of  different  sizes 
and  the  model  assumed  droplets  of  one  distinct  size.  The 
results  reported  in  Table  5*13  show  that  for  the  four  size 
ranges,  each  range  was  retained  in  different  amounts  and 
thus  the  poor  fit  of  the  model  could  be  associated  with  the 
nature  of  the  model.  If  the  emulsion  droplet  size  ranges 
are  retained  in  different  amounts,  then  the  model  should 
allow  for  this  phenomenon.  Also  a  different  form  of  the 
model  to  those  used  in  this  study,  may  produce  a  better  fit 
to  this  type  of  produced  oil  concentration  profile.  A  model 
that  may  produce  a  better  fit  will  be  discussed  in  Section 

5.6.5 

5.6.2  N-Fractional  Size  Emulsion 

The  model  for  a  n-fractional  size  emulsion  is 
described  because  emulsions  are  made  up  of  droplets  of 
different  sizes,  rather  than  of  one  size  and  each  size  has 
a  different  probability  of  capture . 

The  raw  and  processed  data  of  the  produced  emulsion 
droplet  size  distribution  profile  is  reported  in  Appendix 
D.  A  summmary  of  some  of  the  pertinent  data  collected  in 
run  129  is  shown  in  Table  5.13*  The  fit  of  the  model,  for 
each  size  range,  to  the  experimental  data  is  shown  in  Table 
D8  and  Figure  5*21.  The  fit  of  the  model  of  the  combined 
size  range  data,  to  the  produced  emulsion  concentration 
profile,  is  shown  in  Table  D9  and  Figure  5*20.  The  best  fit 
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values  for  the  coefficients  Ho ,  (3  and  e  for  each  size  range 
are  shown  in  Table  5*15* 


TABLE  5.15 

COEFFICIENTS  FOR  N=4  EMULSION  MODEL 


Ho 

13 

e 

Channel 

sec  ^ 

- 

- 

- 

3+4 

0.587 

0.0311 

0.679 

0.0389 

5 

1.162 

0.0244 

0.662 

0.0238 

6 

o.64o 

O.O33 

0 . 666 

0.0251 

7+ 

0.227 

0.048 

0.662 

0.032 

The  results  shown  in  Table  D8  and  Figure  5*21  demonstrate 
that  the  model  can  effectively  describe  the  experimental 
curves  for  the  individual  droplet  size  ranges.  Comparison 
of  the  data  shown  in  Table  D9  and  Figure  "5*20  shows  that 
the  n-fractional  sized  emulsion  model  fitted  the 
experimental  data  better  than  emulsion  flow  model  3*  The 
n-fractional  model  provided  a  better  fit  because  with  the 
nature  of  this  model,  it  could  allow  for  each  size  range  of 
the  emulsion  to  be  retained  in  different  concentrations. 

Examination  of  the  data  shown  in  Table  5*15>  shows 
that  for  all  the  channels,  the  value  of  (3  was  similar  to 
<7m  and  this  further  confirms  the  validity  of  model  3» 
Equation  5- 11 »  over  the  other  two  models  used  in  this 
study.  Although  the  n-fractional  model  effectively  fitted 
the  data  presented  in  the  four  size  ranges,  the  nature  of 
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FIGURE  5.21  EXAMINATION  OF  MODEL  3  FOR  EACH  SIZE  RANGE  :  RUN  129 
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the  kinetic  equation  is  not  of  the  correct  form.  As 
explanied  in  Section  5*5*2,  the  concentration  for  a 
certain  size  range  may  he  greater  during  the  run,  than  when 
steady  state  conditions  are  attained.  As  such,  the  kinetic 
equation  of  the  form 


do 

dt 


f  (  a  ) 


should  be  replaced  by  a  form 

=  H  f ( a )  -  H  f  (  a  ) 

dt  0  r 

where  the  first  term  designates  the  rate  at  which  the 
emulsion  concentration  tends  to  its  maximum  saturation  and 
the  second  term,  the  rate  at  which  the  particles  leave 
their  constriction  sites.  The  combination  of  these  two 
rates,  yielding  the  rate  at  which  the  concentration  tends 
to  its  equilibrium  concentration.  However,  in  these  early 
stages  of  modelling  this  highly  complicated  process,  the 
simplified  model  seems  to  be  to  a  certain  extent 
justified . 

5.6.3  Emulsion  Saturation  and  Concentration  Profiles 

The  plugged  oil  saturation  and  the  flowing  emulsion 
concentration  variation  with  time  and  position,  in  the 
porous  medium,  can  be  evaluated  using  the  emulsion  flow 
model.  Figure  5*22  shows  the  change  in  the  plugged  oil 
saturation  with  time  for  run  129  using  the  n=l  size  ranged 
emulsion  flow  model.  Figure  5*23  shows  the  concentration  of 


124 


^o/d  ‘UOj4U9|0y  SS0|UOISU01UJQ 


FIGURE  5.22  RETAINED  OIL  CONCENTRATION  PROFILE 
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FIGURE  5.23  FLOWING  EMULSION  CONCENTRATION  PROFILE 


126 


the  flowing  emulsion  for  run  129. 

The  properties  exhibited  by  the  model  and  shown  in 
Figures  5*22  and  5-23  are: 

1)  the  retained  oil  concentration  at  the  inlet  end  of 
the  core  can  increase  slowly  with  time  to  its  maximum 
concentration.  This  relationship  is  solely  based  upon  the 
form  of  the  kinetic  equation,  see  Section  3*3.1 . 

2)  the  model  shows  the  continually  changing  shape  of 
the  plugged  and  flowing  concentration  profiles.  If  the 
model  predicts  the  flood  front  is  unstable,  then  the 
concentration  profile  will  continually  change.  If  the 
flood  front  is  stable,  then  the  shape  of  the  concentration 
profiles  will  change  until  steady  state  conditions  are 
achieved . 

A  property  of  the  model  not  shown  in  Figures  5*22  and 
5»23>  is  the  rate  at  which  the  oil  front  moves  through  the 
porous  medium.  For  early  times,  the  model  allows  the  oil 
front  to  move  through  the  system  at  a  rate  greater  than  the 
superficial  velocity.  At  later  times,  when  the  form  of  the 
oil  front  has  developed  further,  the  oil  front  can  fall 
behind  the  position  designated  by  the  superficial  velocity. 
Thus  in  the  laboratory,  the  measured  breakthrough  time  can 
be  a  function  of  the  frontal  development  and  therefore  a 
function  of  the  core  length. 

The  model  should,  however,  only  be  used  to  provide 
concentration  profiles  within  the  core,  when  the  fit  of 
the  model  to  the  experimental  data  is  satisfactory,  because 
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in  the  derivation  of  the  model,  the  material  balance 
equation  was  approximated  to  obtain  a  solution  to  the 
differential  equations.  Thus  when  the  fit  of  the  model  is 
poor,  the  material  balance  for  the  system  is  violated. 

For  further  work  in  this  area,  it  may  be  better  to 
define  the  rate  equation  in  terms  of  the  total  oil 
saturation  in  the  core,  rather  than  the  plugged  saturation. 
In  that  way,  the  material  balance  equation  and  the  plugging 
rate  equation  could  be  solved  without  having  to  apprximate 
the  material  balance  equation.  If  the  model  were  derived  in 
terms  of  the  total  oil  saturation  and  the  flowing  emulsion 
concentration,  then  the  model  could  be  used  to  evaluate: 

1)  the  fractional  flow  curve, 

2)  whether  the  flood  is  stable  or  unstable, 

3)  approximate  time  for  the  flood  front  to  become' 
stabilised  and 

4)  change  of  the  fractional  flow  curve  with  time. 

The  present  models  can  only  provide  an  estimate  of  whether 
the  front  will  stabilise  and  the  time  to  stabilise,  by 
examining  the  relationship  between  arn  and  F  with  time  and 
position 

5.6.4  Pressure  Drop  Model 

The  predicted  pressure  drop  profile  for  run  129  is 
shown  in  Figure  5-24.  The  pressure  drop  across  the  core  was 
predicted  by  numerically  integrating  Darcy's  Law 
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FIGURE  5.24  PRESSURE  DROP  MODEL  FOR  RUN  129 
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(5.10) 


The  relationships  for  the  viscosity  and  permeability  at 
any  position  and  time  were  evaluated  by: 

1)  using  the  model  to  predict  the  flowing  emulsion 
concentration  T  ,  and  the  plugged  oil  concentration 

2)  using  the  viscosity  correlation 

2 

fie  =  /xw  (  l  +  2.666  T  +  5- 333  T  )  and 

3)  using  the  permeability  relationship 

K(  a  )  =  K  -  (  1  -  K  )  o 

w  e  - 

<7 

m 

This  model  was  used  because  each  pore  size  is 
assumed  to  have  an  equal  probability  of  being 
plugged.  Therefore,  the  permeability  reduction 
is  directly  proportional  to  the  plugged  oil 
concentration,  assuming  the  blockage  coefficient 
b  is  constant. 

The  pressure  drop  model  showed  a  satisfactory  correlation 
with  the  experimental  data  until  one  pore  volume  of  emulsion 
had  been  injected.  After  one  pore  volume,  the  pressure  drop 
model  could  not  provide  a  satisfactory  correlation  because 
the  flow  model  on  which  it  is  dependent,  also  provided  a 
poor  correlation. 
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5.6.5  Applicability,  Capability  and  Deficiencies  of  the  Model 

The  emulsion  flow  model  is  obviously  a  rough 
simplification  of  the  whole  set  of  physical  and  mechanical 
phenomena  which  accompany  the  emulsion  plugging  and  flow 
process  in  a  porous  medium.  However,  as  we  are  only  in  the 
initial  stage  of  the  investigation  on  this  highly 
complicated  process,  the  use  of  such  a  simplified  model 
seems  to  be  to  a  certain  extent  justified. 

The  model  can  predict  the  plugged  oil  saturation  and 
the  concentration  of  the  flowing  emulsion  within  the  core 
by  fitting  the  model  to  the  produced  emulsion  concentration 
profile.  The  model  is  only  valid  when  it  provides  a 
satisfactory  fit  to  the  produced  emulsion  concentration 
profile  because  the  material  balance  equation  was 
approximated  to  obtain  a  solution  to  the  differential 
equations.  The  model  overcomes  the  need  for: 

1)  knowledge  of  the  relative  permeability,  in  order 
to  predict  the  concentration  profiles, 

2)  a  capillary  pressure  -  saturation  relationship, 
and 

3)  the  inlet  concentration  changing  from  zero  to  the 
maximum  concentration  at  the  instant  flow  begins. 

This  latter  requirement  is  replaced  by  an  equation  that 
describes  the  saturation  change  in  the  system  is  dependent 
on  the  plugged  oil  saturation  within  the  system  and  the 
concentration  of  the  flowing  emulsion.  This  system  rate 
equation  is  allied  to  the  rate  of  formation  of  a  product 
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in  a  chemical  reaction. 

The  deficiencies  of  the  model  in  its  present  form  are: 

1)  the  model  can  only  predict  the  variation  in  the 
plugged  oil  concentration  and  the  flowing  emulsion 
concentration  with  time  and  position,  and  not  the 
overall  oil  saturation  S  . 


o 

2)  the  correct  form  of  the  kinetic  equation  has  to 

he  found  by  trial  and  error.  Model  3  in  this  study 
provided  a  satisfactory  fit  to  the  early  time  data 
in  runs  119»  120  and  129,  but  not  128. 

3)  the  model  is  only  valid  for  a  linear  system  because 
<7m  is  assumed  constant  along  the  system  length . 

In  a  radial  system,  because  the  pressure  drop  per 
unit  length  is  always  decreasing  with  length  away 
from  the  wellbore,  the  oil  retention  would  be 
increasing  with  length  away  from  the  wellbore . 

4)  the  material  balance  equation  has  to  be 
approximated  to  obtain  a  solution  to  the 
differential  equations. 

The  latter  defiency  could  be  rectified  by  defining  the 
kinetic  rate  equation  in  terms  of  the  overall  oil  saturation 
S  ,  rather  than  the  plugged  oil  saturation.  In  this  way,  the 
material  balance  equation 


£fo 

d  t 


q 

+  - 

A  </> 


dr 

-  =  o 

d  L 


(3-4) 


could  be  used  in  the  model,  rather  than  the  approximated 
form  of  Equation  3-3,  Equation  3*19 
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do  _  dr  q  dF 

-  +  e -  +  -  -  =  o  (3-19) 

dt  dL  A  0  dL 

A  kinetic  rate  equation  that  migth  provide  a  better 
correlation  with  the  experimental  data,  than  model  3  of 
this  study,  is 

2k-  =  H  (  o  -  a  )(  X  -  a  )(  13  +  a  )T  (5-12) 
at  0 

A  comparison  of  the  from  of  the  above  model,  with  model  3 > 
is  shown  in  Figure  5-25*  This  new  model  exhibits  the  same 
high  rate  of  plugging  at  low  concentrations  of  a,  as  does 
model  3  of  this  study.  This  model  might,  however,  provide 
a  better  fit  to  the  experimental  data  than  model  3>  because 
at  concentrations  of  o  approaching  am  ,  the  model  predicts 
a  lower  plugging  rate .  This  lower  plugging  rate  would  increase  the 
increase  the  time  required  for  the  system  to  reach  steady 
state  conditions.  Therefore,  the  new  model  has  the  same 
form  as  model  3>  over  the  range  where  model  3  effectively 
fitted  the  experimental  data,  and  has  an  improved  form 
over  the  range  where  model  3  deviated  away  from  the 
experimental  data. 


Plugging  Rate  da/d t 
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FIGURE  5.25  ILLUSTRATION  OF  THE  FORM 

OF  THE  NEW  KINETIC 


EQUATION 
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CONCLUSIONS 


1 .  Emulsions  can  effectively  reduce  the  permeability  of 
a  core,  with  any  initial  permeability,  if  the 
emulsion  droplets  are  larger  than  the  pore  throats  in 
the  porous  medium. 

2.  There  are  three  degrees  of  core  permeability  reduction, 
zero,  partial  and  complete,  which  are  dependent  on  the 
relationship  between  the  pore  size  in  the  porous  medium 
and  the  droplets  in  the  emulsion. 

3>  When  the  permeability  is  partially  reduced  there  is  no 
simple  relationship  between  the  fractional  permeability 
of  a  core  after  emulsion  flow  and  the  core  properties 
of  porosity  and  permeability.  The  shape  of  the 
fractional  permeability  versus  flow  rate  or  pressure 
drop  curve  is  a  function  of  the  pore  size  distribution 
in  the  porous  medium. 

4.  The  flow  of  an  emulsion  through  a  porous  medium  is  a 
non-Darcy  flow  process  when  the  permeability  of  the 
core  is  reduced  by  the  emulsion  flow.  When  the  emulsion 
does  not  plug  the  core  and  reduce  the  permeability,  the 
flow  is  a  Darcy  flow  process. 

5.  For  a  porous  medium  that  is  plugged  by  an  emulsion,  the 
fractional  permeability  of  the  porous  medium  increases 
with  an  increase  in  the  flow  rate  or  pressure  drop. 

6.  The  parameter  d^/de ,  the  ratio  of  the  sand  grain  size 
to  the  deformed  size  of  the  largest  droplet  in  the 
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emulsion,  can  be  used  to  predict  whether  the  porous 
medium  will  be  plugged  by  emulsion  flow. 

?.  The  capillary  number,  N  can  be  used  to  estimate  the 

c 

flow  rates  over  which  an  emulsion  can  plug  a  porous 
medium . 

8 .  The  concentration  of  the  larger  droplets  retained  in 
the  core  is  larger  than  for  the  smaller  droplets, 
possibly  because  there  are  more  plugging  sites 
available  within  the  core  for  the  larger  droplets. 

9.  The  maximum  concentration  for  a  particular  droplet  size 
does  not  have  to  be  equal  to  the  concentration  at 
steady  state  conditions . 

10.  Model  3  of  the  three  models  investigated  in  this  study, 
provided  the  best  fit  to  the  experimental  data.  This 
model  was  able  to  effectively  describe  the  experimental 
outlet  emulsion  concentration  profiles  for  early  times. 


' 
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RECOMMENDATIONS 


1.  To  evaluate  the  kinetic  equation 

-f^-  =  Ho  (  am  -  a  )(  X  -  a  )(  0  +  a  )T 
where  X  >  am 

/ 3  less  than,  equal  to  or  greater  than  crm 

2.  To  develop  a  kinetic  rate  equation  in  terms  of  the 
overall  oil  saturation,  S  ,  rather  than  the  plugged 
oil  concentration,  o. 

3 .  To  use  a  microwave  attenuation  instrument,  to  obtain 
the  oil  saturation  profiles  within  the  core  during 
emulsion  flow.  In  this  way,  the  model  can  be  compared 
with  both  the  outlet  emulsion  concentration  and  the 
profiles  within  the  core . 

4.  To  use  long  cores,  of  the  order  of  one  metre  long,  so 
that  the  emulsion  front  can  develop  further  and  the 
model  can  be  fitted  to  the  experimental  data  more 
accurately . 
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The  Characteristic  Emulsion  Droplet  Length: 


The  length  of  the  deformed  emulsion  droplet  was 
calculated  using  Equation  5*4 
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The  equation  provides  an  approximate  estimate  of  the  emulsion 
droplet  length  provided  the  droplet  length  is  not  greater 
than  the  diameter  of  the  sand  grain. 

Thus  the  equation  is  probably  valid  for  all  the  runs 
of  this  study  except  runs  100,  103  and  106 .  For  all  these 
runs,  the  deformed  length  is  too  small,  which  leads  to 
high  estimates  for  the  critical  capillary  number,  a  high 
value  for  the  deformed  emulsion  droplet  radius  and  a  low 
estimate  of  the  parameter  d  .  Thus  a  low  estimate  for  the 


length  of  the  deformed  emulsion  droplet,  produces  an 
optimistic  impression  of  whether  a  droplet  will  plug  a  pore 
constriction . 

The  Deformed  Emulsion  Droplet  Radius: 

The  deformed  emulsion  droplet  radius  was  calculated 
using  the  equation 
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which  can  be  written  as 


68942  AP 

x  x  L j  =2x0. 24  x  1 

30.48  AL  a 

in  consistent  units, 

O 

where  68942  dynes/cm  per  psi 
30.48  cm  per  foot. 

Thus  the  equation  solved  by  trial  and  error  was 


4712.25  AP 

Al 


2(  ro  -  re  }  1 

r2  +  r2  +  r  r  J 

o  e  o  e 


where  AP 

-  Pressure  gradient  (  psi  /  ft  ) 

Al 

rQ  Emulsion  droplet  size  (  microns  ) 

rg  Deformed  emulsion  droplet  size  (  microns  ) 

7  Interfacial  tension  (  dynes  /  cm  ) 


Critical  Capillary  Number: 

The  critical  capillary  number  was  calculated  using  the 
equation 

K  2  1  f  1  1  ^ 

[  Nc  J  critical  =  £  73^  —  [-  ~  j 

d  v  p  o  * 

where  K  permeability  (  darcies  ) 

r  pore  throat  size  (  r  /  6.493  ) 
p  ^  g 

rQ  emulsion  droplet  size  (  microns 


(  microns  ) 
) 
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deformed  droplet  length  for  the  system,  calculated 
using  r^  and  rQ  in  equation  (  microns  ) 

porosity  (  -  ) 


APPENDIX  B 


SENSITIVITY  ANALYSIS  ON  THE 


CRITICAL  CAPILLARY  NUMBER 


146 


Sensitivity  Analysis  on  the  Critical  Capillary  Number: 


The  critical  capillary  number  was  calculated  using 
the  equation 

Nc  =  K  2  1  Cl 

4>  1.013  Ld  \rp 

Examination  of  this  equation  shows  that  for  an  emulsion 
droplet  with  a  radius  similar  to  the  pore  throat  radius, 
then  a  small  error  in  either  size  leads  to  a  significant 
error  in  the  capillary  number  because  of  the  sensitivity 
associated  with  evaluating  the  group: 


While,  when  the  emulsion  droplet  size  is  significantly 
different  from  the  pore  throat  size,  a  small  error  in  either 
reading  only  produces  a  small  error  in  the  capillary 
number.  This  phenomenon  is  demonstrated  in  the  analysis 
shown  in  Table  B1 . 

Therefore  when  the  emulsion  droplet  size  and  the  pore 
throat  size  are  similar,  the  capillary  number  is  very 
sensitive  to  the  input  data  and  the  analysis  should  not  be 
used  to  predict  the  flow  rate  at  which  no  plugging  will 
occur.  When  the  droplet  size  and  the  pore  throat  size  are 
significantly  different,  the  errors  associated  with  the 
analysis  are  minimal  and  the  concept  can  be  used  to  predict 
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) 
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the  flow  rate  at  which  no  plugging  will  occur. 


TABLE  B1 

Sensitivity  Analysis  on  the  Critical  Number 


Pore  Droplet  Droplet  Critical  Percent 

radius  radius  length  flow  rate  change 


r 

P 

microns 

r 

0 

microns 

Ld 

microns 

N 

c 

0 

0 

\  & 

^  0 

r 

in  r 

P 

% 

RUN 

104 

12.66 

12.7 

25.44 

8.02 

10 

207 

- 

12.67 

12.7 

25.43 

6.01 

10 

155 

+0.08 

12.68 

12.7 

25.42 

4.01 

10 

103 

+0.16 

12.69 

12 .7 

25-41 

2.00 

10 

52 

+0'.  24 

12.65 

12.7 

25.45 

1.00 

10 

259 

-0.08 

12.6 

12.7 

25.5 

2.01 

10 

518 

-0.47 

RUN 

110 

7.85 

10.5 

23.65 

5.34 

10 

2.1 

10 

- 

7.86 

10.5 

23-64 

5.31 

10 

2.09 

10 

+0.1 

7.84 

10.5 

23.66 

5.36 

10 

2.11 

10 

-0.1 

7.9 

10.5 

23.6 

5.2  2 

10 

2.06 

10 

+0 . 6 

7.8 

10.5 

23.7 

5  .46 

10 

2.15 

10 

-0.6 

Percent 
change 
in  N 

c 

% 


-25 

-50 

-75 

+25 

+150 


-0.2 
+o .  5 

-1.9 

+2.4 
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Cl.l  Outlet  End  Emulsion  Concentration 


The  outlet  end  emulsion  concentration  profile  and 
the  maximum  oil  saturation  for  the  core  were  evaluated  from 
the  emulsion  viscosity  and  the  viscosity  correlation.  Use 
of  the  viscosity  correlation  for  the  series  of  collected 
emulsion  samples  yielded  a  stepwise  outlet  end  emulsion 
concentration  profile  from  which  a  smoothed  line  was  drawn 
to  produce  a  continuous  profile. 

The  maximum  oil  saturation  within  the  core  was 
evaluated  using  the  stepwise  emulsion  concentration  of  the 
produced  emulsion.  This  saturation  was  calculated  by 
applying  a  correction  for  the  volume  of  emulsion  in  the  end 
blocks  to  the  volume  of  oil  retained  within  the  core 
holder . 

The  raw  and  processed  data  for  runs  119»  120,  128  and 
129  is  shown  in  tables  C.l  -  C.3>  C.4  -  C.6,  C.7  -  C.9  and 
C.10  -  C.12,  respectively. 
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RUN  119 


Core  pore  volume  159*6  cc 
Inlet  end  block  volume  6. 3  cc 
Outlet  end  block  volume  7,4  cc 
Outlet  tube  volume  2.3  cc 
Total  fittings  volume  16.0  cc 
Oil  volume  in  the  core  holder  29.7  cc 
Oil  volume  in  the  core  27.3  cc 
Maximum  oil  saturation  in  the  core  0.171 


After  emulsion  breakthrough,  the  pore  volumes  injected 
into  the  core  is  given  by: 

Pore  volumes  injected  Q  -  16.0 

into  the  core  ~  ^ 

where  Q  =  Cumulative  emulsion  injected  into  the  core 
holder  (  cc  ) 

The  stepwise  and  smoothed  outlet  end  emulsion  profile 
is  shown  in  Figure  Cl.  The  experimentally  recorded  data  for 
this  run  is  shown  in  Table  Cl,  with  the  smoothed  outlet  end 
profile  shown  in  Table  C2.  The  theoretically  determined 
profile  evaluated  using  the  model  is  recorded  in  Table  C2. 

A  summary  of  the  pertinent  experimental  data  and  the  values 
for  the  coefficients  of  the  model  are  shown  in  Table  C.3* 
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FIGURE  C.1  OUTLET  END  CONCENTRATION  :  RUN  119 
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cc 

104 

152 
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238 

280 

322 

372 

4l6 

458 

500 

546 

586 

642 

688 

732 

780 

824 

868 

906 

950 
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TABLE  Cl 


Experimental  Data 


Incremental 

Volume 

Emulsion 

Viscosity 

Emulsion 

Cone  . 

cc 

2  -1 

mm  s 

104 

1.0 

0.0 

48 

1.002 

0.001 

42 

1 .23 

0.075 

44 

1.475 

0.141 

42 

1  .478 

0.142 

42 

1 .480 

0.143 

50 

1 .481 

0.143 

44 

1 .483 

0.143 

42 

1 .483 

0.143 

42 

1 .485 

0.144 

46 

1 .487 

0.144 

40 

1 .489 

0.145 

56 

1 .49 

0.145 

46 

1.493 

0.146 

44 

1.495 

0.146 

48 

1  .497 

0.147 

44 

1.498 

0.147 

44 

1.50 

0.147 

38 

1  .501 

0.148 

44 

1.504 

0.15 

Oil  in  the 
Core 
Holder 
cc 

15.6 

22.7 

25.9 
26.3 
26  .6 

26.9 

270 

27.6 

27.9 
28.1 
28.3 
28  .6 

28.9 

29.1 

29.2 

29.4 

29.5 

29.6 

29.7 
29.7 
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TABLE  C2 


Produced  Emulsion  Concentration  Profile 


Pore  Volumes  Dimensionless  Predicted  Emulsion 
Injected  Emulsion  Outlet  Outlet  End 

Concentration  Concentration 


r/i]  r/n 


0.213 

0.0 

0.0 

0.526 

0.0 

0.0 

0.683 

0.0 

0.0 

0.840 

0.060 

0.044 

0.902 

0.140 

0.140 

0.965 

O.338 

0.363 

0.996 

0.550 

0.516 

1.027 

0 . 666 

0 .666 

1.090 

0.857 

0.875 

1.153 

0.917 

0.961 

1.310 

0.945 

0.998 

1 . 466 

0.950 

1.000 

1.623 

0.950 

1.000 

1.779 

0.959 

1.000 

1 54 


TABLE  C3 

Summary  of  Experimental  Data  and  Coefficients  for  Model 


Absolute  core  permeability 

8.12 

darcies 

Core  permeability  to  emulsion 

6.71 

darcies 

Fractional  permeability  to  emulsion 

0.82  6 

Maximum  oil  saturation 

0.171 

Plugged  oil  saturation 

0.0835 

Fractional  porosity  open  to  flow 

0.583 

Blockage  coefficient 

4.99 

Breakthrough  time  (  pore  volumes 

) 

0.72  7 

Coefficients  of  the  Mathematical 

Model 

Rate  coefficient  H 

0 

0.812 

-1 

sec 

/ 3  Factor 

0.0764 

746 

sec  . 

a 

35-6 

cm  ^ 
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RUN  120 


Core  pore  volume  157-3  cc 
Inlet  end  block  volume  6.3  cc 
Outlet  end  block  volume  7-^  cc 
Outlet  tube  volume  2.3  cc 
Total  fittings  volume  16.0  cc 
Oil  volume  in  the  core  holder  30.1  cc 
Oil  volume  in  the  core  27-7  cc 
Maximum  oil  saturation  in  the  core  O.176 


After  emulsion  breakthrough,  the  pore  volumes  injected 
into  the  core  is  given  by: 

Pore  volumes  injected  Q  -  16.0 

into  the  core  “  1 „ 

J-3/  -3 

where  Q  =  Cumulative  emulsion  injected  into  the  core 
holder  (  cc  ) 

The  stepwise  and  smoothed  outlet  end  emulsion  profile 
is  shown  in  Figure  C2.  The  experimentally  recorded  data  for 
this  run  is  shown  in  Table  C4,  with  the  smoothed  outlet  end 
profile  shown  in  Table  C5 •  The  theoretically  determined 
profile  evaluated  using  the  model  is  recorded  in  Table  C5  - 
A  summary  of  the  pertinent  experimental  data  and  the  values 
for  the  coefficients  of  the  model  are  shown  in  Table  C6 . 
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FIGURE  C.2  OUTLET  END  CONCENTRATION  :  RUN  120 
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110 
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232 

274 

326 
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544 

590 

634 

682 

724 

764 

808 

856 

900 
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TABLE  C4 


Experimental  Data 


Incremental 

Volume 

Emulsion 

Viscosity 

Emulsion 

Cone  . 

cc 

2  -1 
mm  s 

110 

1.0 

0.0 

36 

1.0 

0.0 

44 

1.11 

O.O36 

42 

1  .45 

0.131 

42 

1.48 

0.139 

52 

1  .50 

0.144 

40 

1.503 

0.145 

50 

1.505 

0.145 

44 

1  .506 

0.145 

48 

1.507 

0.146 

42 

1.509 

0.146 

46 

1  .51 

0.146 

44 

1.513 

0.147 

48 

1.515 

0.148 

42 

1.517 

0.148 

40 

1.52 

0.149 

44 

1  .521 

0.149 

48 

1.523 

0.149 

44 

1.525 

0.15 

Oil  in  the 
Core 
Holder 
cc 

16  .5 

21.9 

26.9 

27.7 
28  .2 
28  .5 

28.7 

28.9 

29.1 

29.4 

29.5 

29.7 

29.8 

29.9 
30.0 

30.1 

30.1 

30.1 

30.1 
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TABLE  C 5 


Produced  Emulsion  Concentration  Profile 


Pore  Volumes  Dimensionless  Predicted  Emulsion 
Injected  Emulsion  Outlet  Outlet  End 

Concentration  Concentration 


r/r;  r/r; 


0.216 

0.0 

0.0 

0.534 

0.0 

0.0 

0.693 

0.0 

0.0 

0.852 

0.033 

0.059 

0.915 

0.099 

0.123 

0.979 

0.206 

0.238 

1.011 

0.319 

0.318 

1.043 

0.437 

0.410 

1.106 

0 . 608 

0.608 

1.170 

0.752 

0.775 

1.329 

0.891 

0.962 

1.488 

O.925 

0.995 

1.647 

0.942 

0.999 

1 .805 

0.953 

1.000 

N 
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TABLE  C6 

Summary  of  Experimental  Data  and  Coefficients  for  Model 


Absolute  core  permeability 

9.5  darcies 

Core  permeability  to  emulsion 

7.27  darcies 

Fractional  permeability  to  emulsion 

0.765 

Maximum  oil  saturation 

0.176 

Plugged  oil  saturation 

0.0987 

Fractional  porosity  open  to  flow 

0.515 

Blockage  coefficient 

4.91 

Breakthrough  time  (  pore  volumes  ) 

0.745 

Coefficients  of  the  Mathematical  Model 


Rate  coefficient  H 

0 

0.422  sec  1 

/ 3  Factor 

0.098 

^bt 

753-2  sec 

a 

18 . 2  cm  ^ 
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RUN  128 


Core  pore  volume  168.1  cc 
Inlet  end  block  volume  6.3  cc 
Outlet  end  block  volume  7.4  cc 
Outlet  end  tube  volume  2.3  cc 
Total  fittings  volume  16.0  cc 
Oil  volume  in  the  core  holder  35-7  cc 
Oil  volume  in  the  core  33*3  cc 
Maximum  oil  saturation  in  the  core  0.198 


After  emulsion  breakthrough,  the  pore  volumes  injected 
into  the  core  is  given  by: 

Pore  volumes  injected  Q  -  16.0 

into  the  core  ~  ^ g  ^ 

where  Q  =  Cumulative  emulsion  injected  into  the  core 
holder  (  cc  ) 

The  stepwise  and  smoothed  outlet  end  emulsion  profile 
is  shown  in  Figure  C3 •  The  experimentally  recorded  data  for 
this  run  is  shown  in  Table  C7»  with  the  smoothed  outlet  end 
profile  shown  in  Table  C8 .  The  theoretically  determined 
profile  evaluated  using  the  model  is  recorded  in  Table  C8 . 

A  summary  of  the  pertinent  experimental  data  and  the  values 
for  the  coefficients  of  the  model  are  shown  in  Table  C9» 
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FIGURE  C.3  OUTLET  END  CONCENTRATION  :  RUN  128 
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TABLE  C7 


Experimental  Data 


Incremental 

Volume 

Emulsion 

Viscosity 

Emulsion 

Cone  . 

cc 

2  -1 

mm  s 

100 

1.0 

0.0 

47 

1.007 

0.002 

43 

1.062 

0.020 

46 

1.334 

0.099 

46 

1 .419 

0.120 

41 

1 .438 

0.125 

45 

1 .463 

0.131 

51 

1 .480 

0.135 

41 

1 .496 

0.139 

40 

1.510 

0.143 

47 

1.515 

0.144 

45 

1.520 

0.145 

4o 

1.523 

0.146 

60 

1 .526 

0.146 

50 

1.535 

0.148 

40 

1.542 

0.15 

Oil  in  the 
Core 
Holder 
cc 

15.0 

22.0 

27  .6 

29.9 

31.3 

32.3 

33.2 

33.9 

34.4 
34.7 
35.0 

35.2 

35.4 
35-6 
35-7 
35-7 
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TABLE  C8 


Produced  Emulsion  Concentration  Profile 


Pore  Volumes  Dimensionless  Predicted  Emulsion 
Injected  Emulsion  Outlet  Outlet  End 

Concentration  Concentration 


r/r,  r/n 


0.183 

0.0 

0.0 

0.480 

0.0 

0.0 

0.629 

0.0 

0.0 

0.778 

0.05 

0.07 

0.92  6 

0.165 

0.184 

1.075 

0.420 

0.403 

1.224 

0.720 

O.669 

1.372 

0.785 

0.858 

1.521 

0.813 

0.947 

1.670 

0.840 

0.982 

2.265 

0.910 

1.000 

2.860 

0.950 

1.000 

3.354 

O.970 

1.000 

4.049 

0.990 

1.000 

TABLE  C9 


Summary  of  Experimental  Data  and  Coefficients  for  Model 


Absolute  core  permeability 
Core  permeability  to  emulsion 
Fractional  permeability  to  emulsion 


3 .64  darcies 
1.90  darcies 
0.522 


Maximum  oil  saturation  0.198 

Plugged  oil  saturation  0.152 

Fractional  porosity  open  to  flow  O.3O8 

Blockage  coefficient  4.56 

Breakthrough  time  (  pore  volumes  )  0.638 


Coefficients  of  the  Mathematical  Model 


Rate  coefficient  H 

0 

(3  Factor 


a 


O.O33  sec 
O.I36 
3859  sec 
8.5^  cm 
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RUN  129 


Gore  pore  volume  160.4  cc 
Inlet  end  block  volume  8.4  cc 
Outlet  end  block  volume  8.3  cc 
Outlet  tube  volume  2.6  cc 
Total  fittings  volume  19*3  cc 
Oil  volume  in  the  core  holder  30.0  cc 
Oil  volume  in  the  core  27.1  cc 
Maximum  oil  saturation  in  the  core  0,169 


After  emulsion  breakthrough,  the  pore  volumes  injected 
into  the  core  is  given  by: 

Pore  volumes  injected  _  Q  -  19*3 

into  the  core  ^ 

where  Q  =  Cumulative  emulsion  injected  into  the  core 

holder  (  cc  ) 

The  stepwise  and  smoothed  outlet  end  emulsion  profile 
is  shown  in  Figure  C4.  The  experimentally  recorded  data  for 
this  run  is  shown  in  Table  CIO,  with  the  smoothed  outlet  end 
profile  shown  in  Table  Cll.  The  theoretically  determined 
profile  evaluated  using  the  model  is  recorded  in  Table  Cll. 

A  summary  of  the  pertinent  experimental  data  and  the  values 
for  the  coeeff icients  of  the  model  are  shown  in  Table  C12. 
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FIGURE  C.4  OUTLET  END  CONCENTRATION  :  RUN  129 


lui 

olu: 

cc 

104 

149 

204 

253 

300 

346 

391 

436 

481 

531 

576 

622 

668 

713 

759 

80  5 

849 


16? 


TABLE  CIO 


Experimental  Data 


Incremental 

Volume 

Emulsion 

Viscosity 

Emulsion 

Cone  . 

cc 

2  -1 
mm  s 

104 

1.0 

0.0 

4  5 

1.003 

0.001 

55 

1.196 

0.062 

49 

1 .447 

0.126 

47 

1.516 

0.142 

46 

1.525 

0.144 

45 

1.532 

0.145 

45 

1.535 

0.146 

45 

1 .541 

0.147 

50 

1.543 

0.148 

45 

1 .544 

0.148 

46 

l  *5^5 

0.148 

46 

1 .5^6 

0.148 

4  5 

1.549 

0.149 

46 

1.550 

0.149 

46 

1.553 

0.15 

44 

1.554 

0.15 

Oil  in  the 
Core 
Holder 
cc 

15.6 
22.3 
27.2 
28.3 

28.7 
29.0 
29.2 
29.4 
29.6 

29.7 

29.8 

29.9 

29.9 
30.0 
30.0 
30.0 
30.0 
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TABLE  Cll 


Produced  Emulsion  Concentration  Profile 


Pore  Volumes  Dimensionless  Predicted  Emulsion 
Injected  Emulsion  Outlet  Outlet  End 

Concentration  Concentration 


r/n  r/r. 


0.191 

0.0 

0.0 

0.503 

0.0 

0.0 

0.659 

0.0 

0.0 

0.815 

0.120 

0.206 

0.971 

0.400 

0.400 

1.127 

0.700 

0.631 

1.282 

0.830 

0.815 

1.438 

0.895 

0.919 

1.594 

0.930 

0.967 

1.75 

0.950 

0.987 

2.373 

0.975 

1.000 

2.997 

0.985 

1.000 

3.620 

0.985 

1.000 

4.244 

1.000 

1.000 

* 
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TABLE  C12 

Summary  of  Experimental  Data  and  Coefficients  for  Model 


Absolute  core  permeability 

4.14  darcies 

Core  permeability  to  emulsion 

2.26  darcies 

Fractional  permeability  to  emulsion 

0.545 

Maximum  oil  saturation 

O.169 

Plugged  oil  saturation 

0.120 

Fractional  porosity  open  to  flow 

0.325 

Blockage  coefficient 

5.61 

Breakthrough  time  (  pore  volumes  ) 

0.659 

Coefficients  of  the  Mathematical  Model 


Rate  coefficient  H 

0 

0.0296sec_1 

(3  Factor 

0.116 

"hat 

3807  sec 

OL 

7.31  cm"1 
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APPENDIX  D 


EMULSION  DROPLET  SIZE  DISTRIBUTION 

RUN  129 
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TABLE  D3 


Production  History 
4  Fractional  Size  Emulsion 


Cumulative 

Emulsion 

Produced 

Channel 

Channel 

Channel 

Channel  Cone . 

Emulsion 

cc 

3+4 

5 

6 

7  + 

104 

0.0 

0.0 

0.0 

0.0 

0.0 

149 

1 .2 

0.5 

0.1 

0.25 

0.001 

204 

52.2 

45.9 

27.0 

24.1 

0.062 

253 

99-0 

97.0 

73.0 

43.2 

0.126 

300 

104.0 

99.6 

93.0 

66 . 0 

0.142 

346 

99.0 

98  .2 

96.0 

85.0 

0.144 

391 

97.5 

99.6 

99.2 

88.0 

0.145 

436 

- 

- 

- 

- 

0.146 

481 

101 

101 

98  .6 

86.9 

0.147 

531 

99 

100 

102 

91 

0.148 

576 

100 

97 

99 

94 

0.148 

622 

100 

97 

103 

93 

0.148 

668 

99 

99 

102 

94 

0.148 

713 

99 

100 

102 

96 

0.149 

759 

8  05 

Emulsion  Droplet  Size 

Not  Evaluated 
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TABLE  D4 

Inlet  Emulsion  Properties 


Inlet  Emulsion  Concentration 

0.15 

Channel 

3+4  Inlet  Concentration 

0.0603 

Channel 

5  Inlet  Concentration 

0.03525 

Channel 

6  Inlet  Concentration 

0.0294 

Channel 

7+  Inlet  Concentration 

0.02505 

Channel 

Channel 

Size  Distribution  (  microns  ) 

3+4 

0 

3-17 

Channel 

5 

3-17  - 

4.00 

Channel 

6 

4.00  - 

5.04 

Channel 

7  + 

5.04  - 

32.0 

dg/d^  parameter,  each  channel 

Channel  3  +  4 

(TO  - 

20.4 

Channel 

5 

20.4  - 

16.2 

Channel 

6 

16.2  - 

12.8 

Channel 

7  + 

12.8  - 

2.02 
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TABLE  D5 

Retained  Oil  Saturation,  Each  Channel 
Used  in  the  Mathematical  Model 
Smoothed  Data 


Channel  3 


Produced  Fluid 

Volume 

cc 

Incremental 

Volume 

cc 

Emulsion 

Cone  . 

Channel 

Cone  . 

Retained 

Oil 

cc 

104 

104 

0.0 

0.0 

2.9 

149 

4  5 

0.001 

0.013 

4.2 

204 

55 

0.062 

0.505 

4.95 

253 

47 

0.126 

0.99 

4.95 

300 

47 

0.142 

1.0 

4.95 

Oil  Volume  in  the 

Core  Holder 

4.95  cc 

Oil  Volume  in  the 

Core 

4.4  cc 

Retained  Oil  Saturation  in  the 

Core 

0.0275 

176 


TABLE  D5  CONTINUED 

Channel  4 

Produced  Fluid 

Volume 

cc 

Incremental 

Volume 

cc 

Emulsion 

Cone  . 

Channel 

Cone  . 

Retained 

Oil 

cc 

104 

104 

0.0 

0.0 

3-3 

149 

4  5 

0.001 

0.01 

4.8 

204 

55 

0.062 

0.537 

5.6 

253 

47 

0.126 

0.99 

5.6 

300 

47 

0.142 

1.0 

5.6 

Oil  Volume  in  the  Core  Holder 
Oil  Volume  in  the  Core 


Retained  Oil  Saturation  in  the  Core 


5.6  cc 
5.0  cc 
O.O31 
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TABLE  D 5  CONTINUED 


Channel  5 


Produced  Fluid 
Volume 
cc 


Incremental 

Emulsion 

Channel 

Retained 

Volume 

Cone  . 

Cone  . 

Oil 

cc 

cc 

104 

104 

0.0 

0.0 

3.7 

149 

4  5 

0.001 

0.005 

5.2 

204 

55 

0.062 

0.459 

6.3 

253 

47 

0.126 

0.97 

6.3 

300 

4  7 

0.142 

1.0 

6.3 

Volume  of 

Oil 

in 

the  Core 

Holder 

6.3  cc 

Volume  of 

Oil 

in 

the  Core 

5-7  CC 

Retained  1 

Oil 

Saturation  in  the  Core 

0.0353 
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TABLE  D5  CONTINUED 

Channel  6 

Produced 

Fluid 

Incremental  Emulsion  Channel 

Retained 

Volume 

Volume 

Cone  . 

Cone  . 

Oil 

cc 

cc 

cc 

104 

104 

0.0 

0.0 

3.1 

149 

4  5 

0.001 

0.001 

4.4 

204 

55 

0.062 

0.27 

5.6 

253 

47 

0.126 

0.73 

5.9 

300 

47 

0.142 

0.93 

6 . 0 

346 

46 

0.144 

O.96 

6.1 

391 

4  5 

0.145 

0.99 

6.1 

436 

4  5 

0.146 

0.99 

6 . 1 

481 

45 

0.147 

0.99 

6.1 

531 

50 

0.148 

1.0 

6.1 

Volume  of 

Oil  in 

the  Core 

Holder 

6 . 1 

cc 

Volume  of 

Oil  in 

the  Core 

5*5 

cc 

Retained  Oil  Saturation  in  the  Core  0.03^7 
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TABLE  D5  CONTINUED 


Channel  7  + 


Produced 

Volume 

cc 

Fluid 

Incremental 

Volume 

cc 

Emulsion  i 

Cone  . 

Channel 

Cone  . 

Retained 

Oil 

cc 

104 

104 

0.0 

0.0 

2.6 

149 

45 

0.001 

0.002 

3-7 

204 

55 

0.062 

0.241 

4.8 

253 

47 

0.126 

0.432 

5.4 

300 

47 

0.142 

0.748 

5.7 

346 

46 

0.144 

O.796 

6 . 0 

391 

45 

0.145 

0.814 

6.2 

43  6 

45 

0.146 

0.848 

6.3 

481 

45 

0.147 

0.902 

6.5 

531 

50 

0.148 

0.910 

6 . 6 

576 

45 

0.148 

O.916 

6.7 

622 

46 

0.148 

0.922 

6  .8 

668 

46 

0.148 

0.934 

6 . 8 

713 

45 

0.149 

0.964 

6.9 

759 

46 

0.149 

0.98 

6.9 

80  5 

46 

0.149 

0.99 

6.9 

ON 

oo 

44 

0.15 

1.0 

6.9 

Volume  of 

Oil  in 

the  Core  Holder  6.9 

cc 

Volume  of 

Oil  in 

the  Core 

6 . 4 

cc 

Retained  Oil  Saturation  in  the  Core  0.04 
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TABLE  D6 


Summary  of  Experimental  Parameters  -  Each  Channel 


Inlet  Retained  Oil  Retained  Oil  Plugged  Oil 


Channel 

Cone . 

r, 

Volume 

cc 

Saturation 

Sm 

Saturation 

°m 

3+4 

0.0603 

9.41 

0.0585 

O.O394 

5 

0.03525 

5 . 66 

0.0353 

0.0241 

6 

0.0294 

5.56 

0.034 7 

0.0253 

7  + 

0.02505 

6  .4 

0.04 

0.032 

TABLE  D7 


Summary  of  the  Mathematical  Model  Parameters 


Channel 

Inlet 

Cone . 

Plugged  Oil 

Saturation 

Rate 

Constant 

13 

Constant 

Breakthrough 

Time 

i] 

Ho 

sec 

PV 

3+4 

0.0603 

0.0394 

0.587 

0.0311 

3923 

0.679 

5 

0.03525 

0.0241 

1.162 

0.0244 

3824 

0.662 

6 

0.0294 

0.0253 

0.641 

O.O332 

3848 

0 . 666 

7  + 

0.02505 

0.032 

0.227 

0.0485 

3821 

0.662 

■i 
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TABLE  D8 

Produced  Emulsion  Concentration  Profile 


Channel  3  +  4 


Pore  Volumes 

Injected 

Dimensionless 

Emulsion  Outlet 

Concentration 

r/r, 

Predicted  Emulsion 

Outlet  End 

Concentration 

r/r, 

0.191 

0.0 

0.0 

O.503 

0.0 

0.0 

0.659 

0.0 

0.0 

0.815 

0.105 

0.104 

0.971 

0.520 

0.520 

1.127 

0.910 

0.910 

1.282 

0.990 

0.990 

1 .438 

1.0 

1.0 

1.594 

1.0 

1.0 

1.75 

1.0 

1.0 

2.373 

1.0 

1.0 

2.997 

1.0 

1.0 

3.620 

1.0 

1.0 

4.244 

1.0 

1.0 
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TABLE  D8  CONTINUED 


Channel  5 

Pore  Volumes  Dimensionless  Predicted  Emulsion 
Injected  Emulsion  Outlet  Outlet  End 

Concentration  Concentration 

r/r;  r/r, 


0.191 

0.0 

0.0 

0.503 

0.0 

0.0 

0.659 

0.0 

0.0 

0.815 

0.105 

0.105 

0.971 

0.410 

0.410 

1.127 

0.820 

0.805 

1.282 

0.960 

0.961 

1 .438 

0.990 

0.993 

1.594 

0.997 

0.999 

1.75 

1.0 

1.0 

2.373 

1.0 

1.0 

2.997 

1.0 

1.0 

3.62 

1.0 

1.0 

4.244 

1.0 

1.0 
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TABLE  D8  CONTINUED 


Channel  6 


Pore  Volumes 

Injected 

Dimensionless 

Emulsion  Outlet 

Concentration 

r/r; 

Predicted  Emulsion 

Outlet  End 

Concentration 

r/r, 

0.191 

0.0 

0.0 

0.503 

0.0 

0.0 

0.659 

0.0 

0.0 

0.815 

0.105 

0.109 

0.971 

0.270 

0.248 

1.127 

0.470 

0.470 

1.282 

0.700 

0.704 

1 .438 

0.865 

O.865 

1.594 

0.945 

0.945 

1.75 

0.980 

0.979 

2.373 

1.0 

1.0 

2.997 

1.0 

1.0 

3.620 

1.0 

1.0 

4.244 

1.0 

1.0 
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TABLE  D8  CONTINUED 


Channel  7  + 

Pore  Volumes  Dimensionless  Predicted  Emulsion 
Injected  Emulsion  Outlet  Outlet  End 

Concentration  Concentration 


r/ij  r/r, 


0.191 

0.0 

0.0 

0.503 

0.0 

0.0 

0.659 

0.0 

0.0 

0.815 

0.08 

0.08 

0.971 

0.240 

0.254 

1.127 

0.340 

0.340 

1.282 

0.440 

0.437 

1 .438 

0.540 

0.540 

1.594 

0.640 

0.639 

1.75 

0.73 

0.728 

2.373 

0.845 

0.933 

2.997 

0.910 

0.986 

3.620 

0.940 

0.997 

4.244 

0.965 

0.999 

185 


TABLE  D9 


Produced  Emulsion  Concentration  Profile 
4  Channel  Model 


Pore  Volumes  Dimensionless  Predicted  Emulsion 
Injected  Emulsion  Outlet  Outlet  End 

Concentration  Concentration 


r/r{  r/r, 


0.191 

0.0 

0.0 

0.503 

0.0 

0.0 

0.659 

0.0 

0.0 

0.815 

0.118 

0.119 

0.971 

0 .400 

0.396 

1.127 

0.707 

0.704 

1.282 

0.834 

0.835 

1  .438 

0.8  95 

0.895 

1 .59^ 

0.928 

0.929 

1.75 

0.951 

0.950 

2.373 

0.974 

0.989 

2.997 

0.985 

0.998 

3.62 

0.985 

0.999 

4.244 

0.994 

1.0 
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EVALUATION  OF  THE  CONSTANTS  IN  THE  MODEL 
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E.l  EVALUATION  OF  am  FROM  Sm 


The  value  for  the  maximum  retained  oil  concentration 
(  (Jm  )  in  the  core  had  to  he  estimated,  because  it  cannot  be 
measured  directly  without  displacing  all  the  oil  that  is 
flowing  within  the  core. 

For  emulsion  flow,  it  can  be  assumed  that  the  emulsion 
will  plug  off  the  small  pore  throats  and  thus  the  emulsion 
will  have  to  flow  through  the  large  pore  throats.  The 
relationship  between  the  porosity  open  to  flow  and  the 
permeability  was  evaluated  by  assuming: 

1)  a  simple  pore  model  for  the  porous  medium,  which  is 
shown  in  Table  E.l 

2)  the  permeability  of  the  pore  is  proportional  to 
the  radius  of  the  pore  squared  and 

3)  the  tortuosity  of  the  medium  is  independant  of  the 
pore  plugging. 

The  calculated  relationship  between  the  porosity  open  to 
flow  and  the  permeability  is  shown  in  Table  E.l  and  Figure 
E.l.  This  pore  model  is  however  close  to  the  relationship 


Fractional  Permeability 
Plugged 


=  | Plugged  Porosity |  (E.l) 


and  because  the  true  pore  size  distribution  of  the  model  is 
unknown,  the  relationship  shown  in  Equation  E.l  was  used  to 


evaluate  am . 


. 


Permeability  Plugged  f  Kfp 
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FIGURE  E.1 


RELATIONSHIP  BETWEEN  PLUGGED 
POROSITY  AND  PERMEABILITY 
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From  Equation  E.l,  then 


Fractional  Permeability 
Open  To  Flow 


2 

1  -  j  Plugged  Porosity  j- 


therefore 


e 


and  remembering 

Sm 

then 


1  -  (  1  -  e  )2 

1  -  yfl  -  K^T 

+  e  r 

Sm  -  (  1  -  V1  "  Kf  ) 


TABLE  E.l 

PORE  MODEL 


PORE  VOLUME  NUMBER 


SIZE  %  PORES 

1  5  1000 

2  10  500 

3  15  333 

4  20  250 

5  20  160 

6  15  83 

7  10  4l 

5  16 


PORE 

POROSITY  PERMEABILITY 
PLUGGED  PLUGGED 


0.05 

0.002 

0.15 

0.019 

0.30 

0.076 

0.50 

0.213 

0.70 

0.426 

0.8  5 

0.655 

0.95 

0  .864 

1.0 

1.0 

8 
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E.2  USE  OF  THE  COMPUTER  PROGRAM  COMPLX 


In  section  5*6.1,  it  was  pointed  out  that  there  was 
not  one  unique  set  of  constants  which  produced  the  best  fit 
to  the  data.  There  is  infact  a  narrow  range  of  closely 
interelated  values  of  Hq  and  (3  that  give  the  same  fit. 

In  order  to  find  the  best  values  for  the  constants 
and  to  ensure  that  the  computer  program  did  not  converge  on 
a  false  plateau  of  the  least  squares  error  parameter  in  the 
analysis,  the  program  was  initiated  from  a  series  of 
starting  values  forthe  input  variables  and  a  true 
convergence  was  assumed  if  the  program  COMPLX  converged  to 
the  same  values  for  each  variable  from  a  number  of  starting 


locations . 
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